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We study phases and phase transitions that can take place in the newly discovered banana (bow- 
shaped or bent-core) liquid crystal molecules. We show that to completely characterize phases 
exhibited by such bent-core molecules a third-rank tensor T ljk order parameter is necessary in addi- 
tion to the vector and the nematic (second-rank) tensor order parameters. We present an exhaustive 
list of possible liquid phases, characterizing them by their space-symmetry group and order param- 
eters, and catalog the universality classes of the corresponding phase transitions that we expect to 
take place in such bent-core molecular liquid crystals. In addition to the conventional liquid-crystal 
phases such as the nematic phase, we predict the existence of novel liquid phases, including the spon- 
taneously chiral nematic (Nt + 2)* and chiral polar (Vr + 2)* phases, the orientationally-ordered but 
optically isotropic tetrahedratic T phase, and a novel nematic Nt phase with Did symmetry that 
is neither uniaxial nor biaxial. Interestingly, the Isotropic- Tetrahedratic transition is continuous 
in mean-field theory, but is likely driven first-order by thermal fluctuations. We conclude with a 
discussion of smectic analogs of these phases and their experimental signatures. 



I. INTRODUCTION 

Liquid crystals are extraordinary systems in that they 
continue to have a revolutionary technological impact 
and to consistently pose new theoretical challenges of 
fundamental interest. They exhibit a rich variety of 
phases with symmetries intermediate between those of a 
the highest symmetry homogeneous isotropic liquid and 
the lowest-symmetry three-dimensional periodic crystal. 
In contrast to their magnetic and ferroelectric solid state 
analogs, whose ordering is driven by energy-entropy com- 
petition, liquid-crystal phase transitions are of predomi- 
nantly entropic origin. Not unrelated to this is the fact 
that, with one exception of the chiral smectic-C* phase 
0, commonly observed liquid crystal phases are non- 
polar. It is, therefore, not surprising that a recent ex- 
perimental discovery by the Tokyo Tech group Q of fer- 
roelectricity in the liquid crystal phase of achiral bend- 
core (banana) molecules has captured the attention of the 
liquid crystal community. Subsequent light microscopy 
studies by the Boulder group || elucidated the molecular 
organization of the newly discovered phase. They con- 
vincingly demonstrated that (what came to be known as) 
the 2?2 phase jl) is an antiferroelectric smectic-C phase in 
which layers spontaneously break chiral symmetry (with 
chirality alternating from layer to layer) and exhibit polar 
order in achiral molecules. Eight distinct phases of bent- 
core molecules, tentatively labeled B\ to B% have been 
identified |^,^] , though most have not been fully charac- 
terized. Two of the most well studied the B^ and B-j, are 
smectic phases consisting of stacks of fluid layers with 
some internal tilt order, and are special because they can 



be switched with an electric field. A material composed 
of achiral nematogens having a ground state that is ferro- 
electric and homogeneously chiral has also recently been 
discovered. This experimental discovery opens up a 
vast new class of achiral molecules that nevertheless ex- 
hibit ferroelectricity and are therefore of interest to the 
liquid crystal display technology. 

While the study of banana liquid crystals has seen 
substantial experimental strides p]-||||, there has been 
relatively little basic theoretical work on this fascinat- 
ing new class of materials. Brand et al. B presented 
an exhaustive, model-independent classification of the 
symmetry-allowed smectic phases, and Roy et al. pi 
introduced a phenomenological Landau model that pro- 
duces many of the banana smectic phases. There are also 
a number of numerical simulations [fll]-[l5| on systems of 
model bent-core molecules that produce nematic phases 
as well as some of the possible smectic phases. Because 
so far, experimental examples of orientationally-ordered 
but spatially homogeneous phases liquid phases |ll| are 
rare, most of the efforts have focussed on the smectic 
phases of bent-core molecules. Here, we will instead fo- 
cus on spatially homogeneous phases, which we will refer 
to as liquid phases, the understanding and classifying of 
whose phase behavior is in many ways a prerequisite to 
the study of more ordered (e.g., smectic) phases, which 
in addition break translational symmetry. We will take 
advantage of the formal developments and analysis pre- 
sented in this paper for the liquid phases in our studies 
of smectic phases, which we defer to a future publication. 

@ 

The first primary conclusion of our work, which forms 
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the starting point of all further analysis presented here, 
is that a third-rank traceless symmetric tensor order pa- 
rameter T y7c , in addition to the usual nematic Q 4 - 7 and 
vector p l order parameters, is necessary in order to cap- 
ture the orientational order observed in experiments on 
banana molecules. Without introducing such an angu- 
lar momentum L = 3 order parameter, only structures 
which have at least C2 and mirror symmetry, such as e.g., 
the biaxial nematic can be captured, thereby precluding 
a first-principles order parameter description of, for ex- 
ample, the most interesting spontaneously-ordered chiral 
phases. 

As we will demonstrate in great detail, once this 
higher-order order parameter is introduced, a com- 
plex web (displayed in Fig. 0) of possible liquid phases 
emerges and, associated with them, a very rich phase be- 
havior. Many of these phases exhibit exotic symmetries 
summarized in Table 1, including Y>2<i-, D2, and C2 sym- 
metries pH , which have not to our knowledge been pre- 
viously identified in spatially uniform (i.e., liquid) states. 
These new anisotropic liquid phases are distinguished 
by the nature of their T* J ordering. The diversity in 
phase diagram topologies originates from a large number 
of symmetry-allowed transition sequences between many 
of the phases that exhibit some nontrivial combination 
of the p l , and/or T y order parameters. Some of the 
novel orientationally-ordered liquid states that we pre- 
dict are the spontaneously chiral nematic (Nt + 2)* and 
chiral polar (Vr + 2)* phases, an optically isotropic tetra- 
hedratic T phase, and a novel nematic Nt phase, with 
D2d symmetry, that is neither uniaxial nor biaxial, but 
rather exhibits a 4-fold improper (£4) rotational symme- 
try about its nematic axis. 

The paper is organized as follows. In Sec. [n|, we present 
a model of a banana-shaped liquid-crystal molecule. By 
considering mass moments of molecules with this shape, 
we are naturally led to introduce the three important or- 
der parameters, p 1 , Q y and T y7c , that are necessary to 
fully describe anisotropic liquid states into which such 
molecules can macroscopically order. In Sec. [II we 
then catalog all thermodynamically distinct liquid phases 
characterizable by these three order parameters. We or- 
ganize these phases according to symmetry groups under 
which they are invariant and present an exhaustive list 
of phase transition sequences allowed by symmetry. We 
construct a Landau theory of the three order parame- 
ters in Sec. IV and analyze the nature of the complicated 
web of phase transitions that it predicts in Sec. [v|, finding 
full consistency with our general group-theoretic analy- 
sis. In Sec. VI, we briefly discuss possible new smectic 



phases that could result when smectic ordering develops 
in the va riou s liquid phases we identify. We conclude 
with Sec. VII by summarizing our results and discussing 



II. SIMPLE MODEL OF BANANA LIQUID 
CRYSTALS: ORDER PARAMETERS 



As is clear from a chemically and geometrically accu- 
rate schematic of a bent-core molecule, shown in Fig{l], 
the most notable characteristic of banana molecules is 
their "V" shape with bent, (on average) planar and there- 
fore achiral cores. This shape earned such molecules a 
name "bow" -shaped. The molecule is characterized by 
G 2v symmetry, defined by a non-polar direction 1/3 (the 
"string" of the bow) , pointing from one endpoint of the 
"V" to the other and an orthogonal polar axis V\ (the 
bow's "arrow"), pointing to the vertex of the "V", as 
illustrated in Fig. |[ 

We can, therefore, expect the molecule to be charac- 
terized by both even- and odd-rank symmetric, traceless 
tensors with the preferred axis of odd-rank tensors along 
the molecular i>i-axis. We can capture these molecular 
features by a simple three-atom rigid bond model of the 
banana molecule illustrated in Fig. ||. 





their relevance to future studies of smectic phases and to 
experiments. 



NOBOW: X - 152 a C-B4- 155 °C - B2 - 173 D C - / 

FIG. 1. Chemically accurate model of a bent-core (banana) 
NOBOW molecule studied in Refs. [2] and [3] that displays 
Isotropic, B2, B4 and Crystal (X) phases. Instantaneously, a 
molecule can be found in a chiral nonplanar configuration (as 
shown on the right), but fluctuates equally between positive 
and negative chiralities, on average is planar and therefore 
achiral. 

As just discussed, associated with each banana 
molecule a is a body-fixed orthonormal coordinate sys- 
tem with unit vectors (v a ,ii v a,2i v a,z)- Molecular C 2v 
symmetry implies invariance under the reflection opera- 
tions f a .3 — > —1^0,3 and v a .2 — > — and the 7r-rotation 
(about v at \) operation u a ^ — * — ^0,2, ^,3 — > —1^0,3, but 
not under the reflection u a ,i — * —v a ,i- In this body-fixed 
frame, the three atoms' coordinates are given by 
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R a ,i = (acos/?)i/ a ,i , 
Ra,2 = -{asm(3)v a< 3 , 
R tt ,3 = (asin/3)^ Qi3 , 



(2.1a) 
(2.1b) 
(2.1c) 



where 2/3 w 120° is the opening angle of the "V" |15|, as 
shown in Fig. 0, and with the origin located on the u a ^ 
axis, half-way between m masses. As is the case for stan- 
dard nematogens where transitions are driven by entropic 
interactions, we expect that the dominant ordering mech- 
anisms of banana liquid crystals will be associated with 
the shape of the molecule and not with electric dipoles. 
We, therefore, focus on the mass-moment tensors as the 
important order parameters for this problem. That is, 
throughout the paper we will assume that the liquid crys- 
tal ordering is driven by steric interactions and, therefore, 
that it is the mass-moment tensors, rather than charge 
moments, that are the primary critical order parameters. 



-y 





(2.4) 



and a third-rank symmetric, traceless tensor 
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(2.5) 



- HA sl3 <» + 5ih <» + 5hi <»)] ■ ( 2 - 6 ) 

These mass-moment tensors can be expanded in terms of 
complete sets of tensors of the appropriate rank formed 
from the vectors v a ,i, ^a,2 and u a3 : 

Ct, a = ci<i (2.7a) 
Ct = c 23 Ql 3 + c 22 (Ql x - Q% 2 ) (2.7b) 
C^ = c 31 T^ + c 32 rjl (2.7c) 



where 



Q%a = <.y a .a--/ 3 a = 1,2, 3, (2.8) 



Ti j ,t = <i<i<i ~ \(S ij <x + 6 jk <i + S H <i) (2.9a) 

ijk _ 



T, 



^,3^,1^,3 



^,1^,3^,3 



(2.9b) 



FIG. 2. A simple three-atom model of a banana molecule 
and a body-fixed orthonormal coordinate system, capturing 
the molecule's characteristic, achiral Ci v symmetry. 

The lowest order mass moment is just the center of 
mass given, in terms of the body-fixed coordinate sys- 
tem, by 
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V m^Ra 



2m + mi 



mi 



2m + mi 



acos/3 I v a i 



(2.2a) 
(2.2b) 



It is natural to define mass moments relative to the cen- 
ter of mass coordinate R" m . Positions of atom /i relative 
to the center of mass are then v a ^ — R a , M — Rem- The 
second mass- moment tensor relative to the center of mass 
can be decomposed into a scalar (proportional to S lJ ) and 
a symmetric, traceless tensor 
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(2.3) 



The third mass-moment tensor can be decomposed into 
a vector part, 



and 



ci 



2mmia 3 cos f3(— mi + 2m cos 2/3) 



(2m + mi) 2 



c 2 3 = 2ma sin (3 
mmi 



mmi 



C22 



C31 



C32 = 



2m + mi 

a 2 cos 2 (3 , 
2m + mi 

2mm 1 (m 2 + 4m 2 ) 3 3 



a 2 cos 2 (3 , 



(2m + mi) 3 
4m 2 mi 



-a 6 cosT (3 



2m + mi 



a sin (3 cos (3 . 



(2.10a) 
(2.10b) 
(2.10c) 

(2.10d) 
(2.10e) 



There are only two independent symmetric-traceless 
molecular parameters in the set Q l J a because the com- 
pleteness relation, 



(2.11) 



a=l 

implies the constraint 
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(2.12) 



a=l 
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For arbitrary orientation of the axes v a ,a, & given 
molecule will in general exhibit 5 independent symmetric- 
traceless tensor mass-moments. However, because we 
have chosen these axes to be symmetry axes of the 
molecule, there are only two independent tensors. Sim- 
ilarly, a given molecule will in general be characterized 
by 7 independent third-rank mass-moment tensors. By 
utilizing 3 rotational degrees of freedom (e.g., Euler's an- 
gles) it is always possible to choose the axes v a , a so that 
there are only four independent such tensors. Our model 
of a bent-core molecule is sufficiently simple, that, with 
our convenient choice of basis vectors i> a ^ a , each molecule 
is characterized by only 2 non-vanishing third-rank ten- 
sors. 

The potential energy of interacting bent-core molecules 
can be expressed in terms of the generalized tensors v l a 1 , 
Q l J a and T£> k and higher-rank tensors. In the phe- 
nomenological treatment we will pursue, it is convenient 
to introduce coarse-grained field versions of these ten- 
sors: 



a 

T^(x) = i^7^(x-x Q ). 



(2.13a) 
(2.13b) 
(2.13c) 



where x Q is the position in the lab frame of the center of 
mass of molecule a and p is the molecular number den- 
sity. Thus, a theory for our model bent-core molecules 
that includes all tensor order parameters up those of 
third-rank would include a single vector order param- 
eter derived from the third-rank mass- moment tensor, 
two second-rank tensors and two third-rank tensors. To 
simplify our discussion, we will consider phenomenologi- 
cal theories with only one second-rank tensor, which we 
denote Q lJ , and one third-rank tensor, which we denote 
T ij . This, however, is not a restriction on our model, 
because if the original theory had all four second- and 
third-rank tensors, we could for example integrate out 
Q % J and order parameters to obtain our model as 
an effective theory, depending only on p 1 , Q v = , and 
rpijk = rpij or( j er parameters. 

Each of these tensors can be expressed in terms of its 
components relative to a space-fixed orthonormal basis 
(111,112,113) = (m, 1, n), with m x 1 = n. To this end, 
we introduce second- and third-rank symmetric-traceless 
orthonormal basis tensors, Jff and Pj k , that transform, 
respectively, under L = 2 and L = 3 representations of 
the rotation group in three dimensions: 
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3/2(7 

\j2(m l m> - PF) , 
Tfl{n l m? ; +mV) 



(2.14a) 

(2.14b) 
(2.14c) 



Jl 3 = y/lfifa'P +fn j ) , 
jj? = y/l/2{mH j + tm?) 



(2.14d) 
(2.14e) 



and 
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1 

T ijk 
'2 

T ijk 
'3 

T ijk 



<r 5/2[n i n j n k - -{5 ij n k + V k n l + 5 kl rt)] (2.15a) 
5 

(m l m j m k - m i l j l k - m?l k t - m k Pl j ) (2.15b) 

(lHH k - l l m?m k - Pm k m l - Z fc mW) (2.15c) 



y/5/12 [m i v?n k + m?n k n l + m k n L n? 

- -(rrfS^ +m j 5 ik +m k 5 ij )] 
5 

4 jk = y^5/12[l i n j n k + Vn k n l + Z fe nV 

_ ~Q i 8 jk +l j 6 ik + l k 6 ij )] 
5 

jijk = _L{ n *( m j m k _ + n i( m i m k _ 
V6 



n k {m l m? - fl j )} 
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- : 7 =(n i m j l k + n l Vm k + m l l J n k + rrinH 



T ijk _ 

+ l i n j m k + l i m :i n k ). 
These tensors are normalized so that 

V p3 k r jk = 5 , 



(2.15d) 
(2.15e) 

(2.15f) 

(2.15g) 

(2.16a) 
(2.16b) 



We can now express our order parameters fields, Eq. 2.13 
in terms of these bases 



rpijk _ TpPik 



(2.17a) 
(2.17b) 
(2.17c) 



The space-fixed orthonormal basis (m, 1, n) can be ro- 
tated to diagonalize the tensors p l , and T y . In 
general, there is no reason why the rotated bases of 
these three order parameters should coincide. We should, 
therefore, in general introduce three bases (111^,14,114) 
where A 6 {p, Q, T}. Any of these bases are fully speci- 
fied by three angles, and we can choose them to eliminate 
up to three components of the tensors Q y and T lJ . In 
particular, we can choose the two independent angles in 
wq to eliminate Q3 and Q4. This leads to 



Q» = S(n* Q n 3 Q --^) 

+B x {m i Q m? Q - l l Q l 3 Q ) + B 2 {w} Q V Q + 1} 



(2.18) 
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The independent angle defining the direction of the pair 
(hiq, \q) can be used to eliminate either B\ or B 2 . A 
similar line of arguments allows us to choose the basis 
(m<r, It, n-r) so that T3 = T4 = T5 = 0. We can parame- 
terize the four remaining components of T^ k in terms of 
an amplitude T and 3 angles, 9i, 62, and 63 and write 



T ijk = T(coa9ili jk + 
+ sin 6*i sin 9 2 



sin 9\ cos 02l2 k 

sin Q\ sin 02 sin O^ltf 



(2.19) 



in the basis (mj, I7 1 , ny). There are other representa- 
tions of the general tensor T 1 - 7 involving other sets of 
four of the tensors I?? k and a different set of three an- 
gles. The representation of Eq. (2.19) is the most useful 
for out purposes. Finally, we can choose the p-basis so 
that 



P = pn p . 



(2.20) 



In what follows, we will, unless otherwise specified, ex- 
press all quantities in the basis (hiq, 1q, uq) = (m, 1, n) 
that diagonalizes with B2 — 0. We will then have to 
worry about the possibility of 7 independent components 
of T^ k and three independent components of p l in this 
basis rather than the angles of the T and p bases rela- 
tive to those of the Q-basis. When Q lJ is zero, we can of 
course choose the p or the T basis. 



III. PHASES AND THEIR SYMMETRIES 



We have just seen that phases of banana molecules can 
be characterized by vector and second- and third-rank- 
tensor order parameters. Before developing a Landau 
field theory for these order parameters and analyzing it 
in mean-field theory, we summarize in this section the 
various phases and their symmetries that can arise from 
these order parameters, and we review possible phase se- 
quences produced by the simplest version of the Landau 
theory. 

Before cataloging the possible phases of our model and 
their symmetries, we observe that many of these phases 
can be successfully described in terms of effective theories 
that are functions of only one of the order parameters 
p i , Q u , or T^ k . Some of the phases, however, partic- 
ularly those of lower symmetry, require two or more of 
these order parameters for their full description. Further- 
more, these order parameters are coupled by rotationally 
invariant contributions to the free energy like p l p J Q VJ , 
QiiT ikl Ti kl , or //(/•" / '-\ and the existence of one kind 
of order induces another. Thus, for example, a model 
based on T^ k alone would miss the fact that the lower- 
rank tensor Q y is automatically induced by T 2 order. 



Phase 


Symmetry 


Order parameters 


V 


Coot, 


P3, S, Ti 


N 




S 


N + 2 


D 2 h 


S, B1.2 


N + 3 


D 3 h 


S, T 2j 3 


T 


T d 


T 6,7 


N T 


D 2(i 


■S) ^6,7 


(N T + 2)* 


D 2 




V + 2 




Pq, S. Bi , Ti , Tr\ 

x yJ 1 1 J-7 i-i U; 

or pi, 5, Bi, T 2 , Ta 


V + 3 




77q S Ti To q 


(V T + 2)* 


c 2 


P3, S, Bi, Ti, Tq, T 7 ; 
or pi, S, Bi, T 2 , T 4 , T 5 


N + V 




Pi, P3, S, Bi, Q 3 , 

Tl, T2, T4, Tq 



Table 1: Anisotropic liquid phases of banana-shaped 
molecules, their symmetries in the Schoenflies notation, 
and their nonvanishing order parameters. The notation 
B\^2, etc. is explained in the text. Some of the phases, 
such as the N + V phase can be characterized by other 
sets of symmetry equivalent order parameters, involving, 
for example, linear combinations of pi andp2 rather than 
Pi alone. 

Table 1 lists the phases we consider, their symmetries, 
and the nonvanishing order parameters that character- 
ize them. This list includes phases with all symmetries 
that can be constructed from the order parameters p % , 
, and T^ k except for the lowest-symmetry phase with 
C\ symmetry ]33||34j], which we do not consider. All 
other point-group symmetries including cubic, icosahe- 
dral, simple tetrahcdral (T), and even lower symmetries 
such as S2, C3, and C2h cannot be characterized without 
the introduction of 4th- or higher-rank tensor order pa- 
rameters. As is customary, we denote the Isotropic phase 
by / and the Nematic phase by N. The N phase has D^h 
symmetry, and it is completely characterized within the 
space of p 1 , Q* 3 ', and T % i k by the single uniaxial order 
parameter S. In general, the N phase will also have non- 
vanishing components of all even rank tensors (explicitly 
induced through Tr(Q"T2 n ) coupling), but we will ig- 
nore them, focusing on the nontrivial order parameters 
of rank 3 or less that actually drive the ordering transi- 
tions. There is a phase with vector or, equivalently, Coot, 
symmetry, which we denote by V. The predominant or- 
der parameter of this phase is the vector p, which we 
take to be along n (i.e., nonzero p 3 ). Once p 3 orders, 
it explicitly induces S and T\ order parameters, through 
the p % p>p k T % i k and p l pr l Q % i couplings, respectively. 

There are a number of phases in which anisotropy de- 
velops in the plane perpendicular to n or p. As we 
will find in Sec. [v], phases that break uniaxial symme- 
try (isotropy of the transverse plane), will exhibit 0(2) 
invariance corresponding to internal rotation within pairs 
of order parameters. For convenience we will collectively 
refer to these pairs as: pi, 2 = (pi,P2), Bi,2 = {Bi,B 2 ), 
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T 2i3 = (T 2 ,T 3 ), 7^5 = (T 4 ,T 5 ), T 6i7 = (T 6 ,T 7 ). Also, 
following Toner |L{|, we denote the phases with r-fold 
anisotropy (or equivalently r-atic order) in the plane 
perpendicular to n by N + r and those with similar 
anisotropy perpendicular to the vector axis p by V + r. 
There is the standard biaxial nematic or TV + 2 phase 
with T)2h symmetry and S and order. There is an 
N + 3 phase with D^h symmetry and nonvanishing S and 
72,3 order. The V + 2 (equivalent to the N + 1 phase) 
and V + 3 phases have C 2v and C3,, symmetry, respec- 
tively. The V + 3 phase develops from the N + 3 phase by 
developing vector order along the n axis. It, therefore, 
has P3 and T\ order in addition to the S and T2,3 order 
of the N + 3 phase. The V + 2 phase has one 2-fold axis 
and two perpendicular reflection planes. Which order 
parameters describe this phase depends on whether the 
vector order p lies along n (which diagonalizes Q y ) or 
perpendicular to n, within the m-1 plane. If p is parallel 
to n, the V + 2 phase is characterized by nonvanishing 
P3, S, B±, T\, and Tq or by a symmetry equivalent set 
such as P3, S, B2, Ti, and T7, where it is understood 
here that T7 is zero if B\ is nonzero and Tq is zero if B 2 
is nonzero (otherwise spontaneous chirality develops, as 
discussed below); if p is perpendicular to n, then it is 
characterized by nonvanishing 5*, B\, T 2 , and T 4 or 
symmetry-equivalent order parameters. 

Schematic representations of the N, N + 2, N + 3, V, 
and V + 2 phases derived from bent-core molecules are 
shown in Fig. |^. The distribution of molecular angles in 
the V + 3 phase is difficult to depict in the format of Fig. 
|3|. Because of the symmetry of the bent-core molecule 
under 1/3 — > —1/3, it is impossible to produce V + 3 sym- 
metry if the molecular i>\ axis is rigidly aligned along p. 
To produce such 3-fold symmetry, it is necessary for the 
molecular U\ axis to be tilted away from the the p axis 
and for its projection onto the plane perpendicular to p 
to have three-fold symmetry. 

A comment about how Ti^, order describes iV + 3 (and 
V + 3) symmetry is useful, r-atic order is generally de- 
scribed by an order parameter of the form (e lr ^), where 
(j) is the angle between a molecular axis in the m-1 plane 
and the m axis. To represent order in this way, we 
introduce the circular basis vectors 



N Phases 



1 

72 



(m±il) 



satisfying 



1, 



(3.1) 



0, (3.2) 



and reexpress I%* k and I l J k as 



T ijk _ 1 / jijk . jijk\ 



T ijk 



V2 





V, 



(b) 

N+2 Phases 

Ak 



(c) 




m 



(e) 

N+3 Phases 



(0 





(h) 





V+2 Phases 

i m,p 



M 



0) 



(k) 



FIG. 3. Schematic representation of the N, N+2, N+3, V, 
and V + 2 phases. Three versions, (a), (b), and (c), of the iV 
phases are depicted with respective predominant alignment of 
1^3, vi, and v-i along n, specifying the direction of the princi- 
pal axis of Q lJ with the largest eigenvalue. The iV + 2 phases, 
(d), (e), and (f), are obtained, respectively, from the iV phases 
(a), (b), and (c) by restricting rotations in the plane perpen- 
dicular to n to have two-fold symmetry, whereas the TV + 3 
phases, (g) and (h), are obtained by restricting these rota- 
tions to have a three-fold symmetry. In the uniaxial V phase 
(i), the molecular v\ aligns on average along p||n, sampling 
equally all orientations about the p axis. The V + 2 phase can 
be produced either by introducing biaxial order perpendicu- 
lar to p and n into the V phase (j) or by introducing vector 
order into the N + 2 phase by aligning p along m (k). 



(3.3) 



G 



Then T 2 I l 2 jk + T 3 I^* = T + I% k + T_/!_ J ' fe , where 



T± = —f= (T 2 =F *T 3 ) = 444T« fc . (3.4) 



If 1/3 is aligned along n, then T± = (e ±i3(?i ) where </> is 
the angle between V\ and m. When 1/3 is not aligned 
along n, the situation is similar though more compli- 
cated. Thus, nonzero describes triadic order in the 
plane perpendicular to n. We will also briefly encounter 
even lower symmetry phases, in which, in contrast to the 
N + r and V + r phases discussed above, the additional 
order develops in a plane that is not perpendicular to the 
established nematic or vector axis. One prominent ex- 
ample is a phase in which the nematic and polar orders 
are neither parallel nor perpendicular. We will refer to 
this Ci/j-symmetry phase as N + V, emphasizing its dis- 
tinction from the N + l (= V + 2) phase, discussed above. 
Although N + V phase (and its N + r + V analogs) is 
quite unlikely to develop in the liquid state, such order 
can quite naturally appear in the smectic-C environment, 
where the additional axis is defined by the smectic layer- 
normal N. 




T N T 1 

FIG. 4. (a) A tetrahedron, that exhibits symmetry identi- 
cal to that of the Tetrahedratic T phase, with all three axes 
of the cube, n, m, 1 equivalent. The T phase can be visualized 
as being composed of banana molecule tetrahedral complexes 
on average decorating edges of randomly positioned but ori- 
entationally-ordered tetrahedra as shown in Fig. |H](a). (b) A 
tetrahedron uniaxially distorted along the n axis, exhibiting 
symmetry of the Nt phase, that is distinguished from the T 
phase by the nonzero nematic order parameter S. A depiction 
of this phase in terms of banana molecules is shown in Fig. 
|b) 




(b) (c) 

FIG. 5. Schematic representation in terms of banana 
molecules of (a) the T phase, (b) the Nt phase, and (c) 
the (Nt + 2)* phase. In the T-phase, bent-core molecules 
align locally with their v% axes aligned on average along the 
six edges of a tetrahedron and their i>\ axes aligned parallel 
to the normals ±m, ±1, and ±n to these edges. Opposite 
edges are (say with normals along n and — n) are orthogo- 
nal so that molecules aligned along opposite edges have per- 
pendicular 1^3 axes. The Nt phase is obtained from the T 
phase by a uniaxial distortion along one of the cubic axes 
as shown in Fig. ^ to favor one pair of crossed bent-core 
molecules over the other two orthogonal pairs. Note the in- 
variance of the Nt phase under the 4-fold improper rotation 
£4 : m — > 1, 1 — » — m, n — » n. The chiral nonpolar (Nt + 2)* 
phase is obtained from the Nt phase by rotating the two 
molecular planes away from 90° to an angle < 5 < n/2 to 
remove the S4 symmetry element, as also illustrated in Fig.H. 



There is only one phase in which T yfc has a nonvan- 
ishing component and in which both p l and Q lJ are zero. 
This phase, which we denote by T and call tetrahedratic, 
has tetrahedral symmetry and is invariant under all 24 
operations of the tetrahedral group T^ p8l . It is charac- 
terized by an arbitrary nonvanishing linear combination 
of Tq and T7, i.e., by the Tq 7 order parameter, and is 
illustrated in Figs. [|(a) and||(a). 



A uniaxial distortion along one of the 3 two-fold tetra- 
hedral axes reduces the T^ symmetry of the tetrahedratic 
phase down to T)2d symmetry. We denote the resulting 
nonpolar phase with this symmetry by Nt- It is charac- 
terized by nonvanishing nematic S and Tq j order param- 
eters, but it is neither a uniaxial nor a biaxial nematic, 
with Tq j breaking the isotropy of the plane transverse to 
n, as illustrated in Figs. ||(b), ||(b), and||(a). 
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front: side: 



top: 
front: 



bottom: 




side: 



V 



FIG. 6. Various representations of the achiral Nt (a) and 
the chiral (Nt + 2)* (b) phases. The two phases are distin- 
guished by their opening angle 8. In the achiral Nt phase, 
S = 7r/2, whereas in the chiral (Nt + 2)* phase, < 5 < n/2. 




parameter, which characterizes the Nt phase. More con- 
cretely, for the choice of the basis m-1, such that Nt is 
exclusively described by nonvanishing S and T7 order pa- 
rameters, the polar achiral V + 2 and the nonpolar chiral 
(Nt + 2)* phases emerge when B 2 and B\, respectively, 
order; equivalently, if it is the nonzero S and T@ that are 
exclusively used to describe the Nt phase, then the roles 
of B 2 and B\ are reversed and transitions to V + 2 and 
(Nt + 2)* take place when B\ and B 2 , respectively, be- 
come nonzero. The polar chiral (Vr + 2)* phase emerges 
from the nonpolar chiral (Nt+2)* phase via development 
of polar order p^ along the existing ncmatic n axis. 








T 2 (S) 










T:7 d 




N+V:C 



l/i 



(V T +2)*:C 2 



The V + 2 phase with C 2v symmetry naturally emerges 
from the Nt phase through the development of longitudi- 
nal polar order p = P3I1 along the existing nematic axis. 
As we will see in Sec.0, once P3 develops in the presence 
of Tq^, a biaxial order Bi 2 with principal axes parallel 
to those of T^ri 1 is explicitly induced. Our final two 
phases, which we denote by (Nt + 2)* and (Vt + 2)*, re- 
spectively, have D2 and C2 symmetry. They are unique 
in that they are spontaneously chiral phases. The non- 
polar chiral (Nt + 2)* phase depicted in Figs. ||(c) and 
|6|(b) is formed from the nonpolar achiral Nt phase by 
the development of biaxial Bi t % order (but in contrast 
to the polar achiral V + 2 phase) with principle axes ro- 
tated exactly by 7r/4 relative to those of the T^n 1 order 



FIG. 7. A flowchart of phase transitions between liquid 
crystal phases illustrated in Fig. ^| Order parameters, which 
become nonzero at each of the transitions and their symme- 
try groups are indicated. For transitions that we have studied 
in detail, we have also indicated the secondary (explicitly in- 
duced by nonlinear couplings) order parameters by placing 
them in parenthesis. 

Alternatively, a transition to it can also take place from 
the polar achiral V + 2 phase by spontaneously breaking 
chiral symmetry via development of biaxial B\ i order 
with principle axes rotated exactly by 7r/4 from those of 
Tg"^ n l order parameter. Not surprisingly, in all the polar 
phases the T\ order parameter is also explicitly induced. 
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Since both the (N T + 2)* and V T + 2)* phases are chi- 
ral, their ground-state configurations will exhibit spatial 
modulations like those of cholesteric and blue phases of 
chiral mesogens. 

Given the achirality of the bent-core molecules, the 
transitions from N T to (N T +2)* and V+2 to (Vr+2)* are 
ones in which chiral symmetry is broken spontaneously, 
and they are, therefore, relevant to the physics of chi- 
ral banana phases. || Figure ^ summarizes the phases 
we treat and the symmetry-lowering transitions among 
them that can take place, i.e., allowed by symmetry 
(as opposed to energetic) considerations. The fact that 
C2v is a subgroup of Coot, {C2v C Coo„) implies that 
there can be a V — > V + 2 symmetry-lowering transi- 
tion. The subgroup structure D 2 C D 2( j C C 0(3) 
and C2v C T>3h C 0(3), where 0(3) is the full orthog- 
onal group including inversions in three dimensions, im- 
ply, respectively that the I -> T — >• N T — >• (N T + 2)* 
and I—fN + 3^V + 2 phase sequences are possible. 
Other phase sequences shown in Fig. ^ follow from simi- 
lar group-theoretic arguments, and are supported by the 
detailed analysis of the Landau mean-field theory given 
in Sec.0. 



IV. CONSTRUCTION OF A LANDAU FIELD 
THEORY 

Having identified the possible spatially homogeneous 
but anisotropic phases |20[| of systems described by first, 
second, and third-rank tensor order parameters, we now 
turn to the study of phase transitions among them. To 
this end we begin by constructing a Landau free energy 
that will describe transitions from the Isotropic phase. 
The appropriate Landau free energy functional is a ro- 
tationally invariant power-series expansion in the order 
parameters p i , , and T % i k . The most general Landau 
free-energy density is produced by sums of scalars formed 
from the tensors p 1 , Q v ~ and T Ij7 \ It can be decomposed 
as 

/ = fp + fQ + It + fpQ + fpT + IqT + fpQT + fq 2 X 2 , 

(4.1) 

where f p , /q, fx are, respectively, the Landau energies 
for independent vector, second-rank, and third-rank ten- 
sor order parameters and the other energies are couplings 
between these order parameters. The vector energy f p is 
given by the standard 0(N = 3) model: 



f p = t Kp (d jP l )(d jP l ) + \r p p l tf + u p (p l p 



(4.2) 



The purely 2nd-rank-tensor part of / is the well-known 
Landau energy for Isotropic-Nematic transition pl"| , p2| , 
given by 
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fQ = \K Q {d k Q^d k Q^) + -r Q Q^Q^ 
-WQQ lj Q 3k Q ki + u Q (Q ij Q i: >) 2 , 



and the purely 3rd-rank-tensor part of / is given by 
f T = \K T {diT vjk diT tik ) + \r T T llk T i]k (4.4) 

In the above expressions we have suppressed the position 
dependence of order parameters, have used an Einstein 
convention for the repeated indices, and have left out 
the dipolar- like ("space-spin" coupling) gradient terms, 
diQ^djQ^, and d h T l ^ k d l2 T l ^ k , that couple internal 
indices of Q % i and T* J * to that of the spatial coordinate 
x. Although this last simplification might modify the 
asymptotic nature of the phase transitions, it obviously 
will not affect our mean-field discussions, valid outside 
of a (typically) narrow critical region. The parameters 
r p — T - T p , r Q ~ T - Tq, and r T — T - T T van- 
ish at the temperatures T p , Tq and Tt, which are de- 
termined predominantly by the interaction potential be- 
tween molecules, that characterize mean-field limits of 
metastability. In writing down quartic nonlinearities in 
/q and /t, we have used the non-obvious relations (valid 
in 3 dimensions) 



1 



{QVQ^f = Q lj Q 3k Q kl Q H 

_^rpijkrpijk^2 _ rpijik 1 rplj 1 k 1 rpij 2 k 2 rplj2k2 



+T l 



'»l»2*3J 1 «l»4*5 7 1 «2*4*6T ,I 3 l 5»6 



(4.5a) 



(4.5b) 



to reduce the number of quartic couplings in /q from two 
to one and in fx from three to two. The lowest-order 
contributions to the coupling energies are 



fpQ = 


-w pQ p l p 3 Q %3 


(4.6a) 


fpT 


-W pT 


(4.6b) 


jo 1 = 


-w Q x Q^ 2 T^ k T Mk 


(4.6c) 


fpQT 


-w pQT P l Q 3k T llk 


(4.6d) 


e 

JQ 2 T 2 — 


-Wi Q^ l Q l2l T lljk T t2jk 


(4.6e) 


f (2) 

J Q 2 T 2 ~ 


—W2 Q n ^ 1 Q' l2 ^ 2 T li: ' lk T l2 ^ 2k 


(4.6f) 


A3) 

JQ 2 X 2 — 


—W3 Q lll2 Qi 1 i 2 T n: ' lk T l2: ' 2k 


(4.6g) 



An) 

,1 JQ2 T 2 - 



The 



(4.3) 



where we have decomposed fcpx 2 as J2 r . 
term p % Qi k T % i k deserves special attention. If the product 
Qjkrpijk j g nonzero t ms term will induce vector p % order. 
Thus, it is possible to have a transition from the TV phase 
that appears to be driven by T yfe but which nonetheless 
develops vector order. In other words, a model expressed 
in terms of Q 1 ^ and T yfc , only, would miss the develop- 
ment of vector order, which by itself is unlikely to order 
in a realistic liquid crystal. 

As usual, the average properties are computed by inte- 
grating over order parameter configurations with a Boltz- 
mann weight with an effective Hamiltonian Ti. = f d 3 xf 







(O) = — J VfVQ^VT^ 0(p\Q ij ,T ijk ) e ~ n/kBT , 

(4.7) 

where 0(p l , Q tJ , T tjk ) is a function of the order parame- 
ters p 1 , Q tJ and T t]k and Z is the partition function. 



the transitions from the / phase are, however, first-order, 
and which of the possible transitions actually takes place 
depends on higher-order terms in the free energy. We 
will thus consider each transition separately. 



V. PHASE TRANSITIONS 

In the preceding two sections, we defined order param- 
eters, we identified possible rotationally anisotropic, but 
spatially homogeneous thermodynamic phases of bent- 
core molecules, and we constructed a Landau theory to 
describe phase transitions among these phases. In this 
section, we will show how each of the phase transition 
sequences depicted in Fig. [?] arises in mean-field theory. 
We will organize our discussion by considering sequen- 
tially each symmetry-lowering transition from each of the 
phases in Fig. ^. Thus, we will first discuss transitions 
from the / phase to the V, N, N + 3, and T phases. 
We will then study transitions from each of these phases, 
that is from the V, N, N + 3, and T phases into lower- 
symmetry phases and so on until the lowest-symmetry 

V + 3, (Vt + 2)*, and N + V phases are reached. 

In our discussion of transitions out of various partially 
ordered phases, it will be useful to display explicitly the 
part of the free energy functional / that couples the order 
parameters p l , Q t3 \ and T 13 , which identify our phases. 
This coupling energy part, which we collectively call f c , 
arises from the Q 3 and Q A parts of fq , from the T 4 parts 

of f T , and from f pQ , f pT , Jqt, fpQT, /qV^, Jq2 T 2, and 

(3) 

/q2 T 2- It determines those order parameters that are 
coupled at harmonic order when long-range order has 
been established in a subset of order parameters, and it 
is, therefore, essential in establishing the nature of the 
phases and phase diagrams for our system. 

A. Transitions from the Isotropic phase 

The Isotropic phase is the phase with the highest 
(0(3)) symmetry. A symmetry-lowering transition to the 

V phase with vector symmetry takes place with the de- 
velopment of p order and one to the N phase with the 
development of uniaxial Q lJ order. As Figs. |?] and [l^ 
indicate, the development of T tjk order in the Isotropic 
phase /, can lead to two distinct phases: the tetrahedral 
phase T with Tq^ order and the N + 3 phase with uni- 
axial Q y order in addition to T2,3 order. At long length 
scales, small fluctuations within the Isotropic phase are 
described by a harmonic free energy density 

fW = \{r p p l p l + tqQVQV + r T T^ k T^ k ). (5.1) 

Thus, which of the fields p % ,Q % - \ or T ljk first becomes 
unstable is determined by which of the set of parame- 
ters Si = {r p ,rQ,rT} first passes through zero. Some of 



1. I — > V transition 




FIG. 8. A schematic for the I —> V transition. 

The I — V transition is driven by the development of 
p order. Since p is a vector, this transition is in the 
well-known 0(3) universality class; it can be described 
in terms of an effective theory involving p only. Below 
this second-order transition, we can take p to point along 
the 3 direction with 

P3 ~ \r P - r pc f°w , (5.2) 

with 0o(3) ~ 0.366, where r pc is the value of r p at the 
critical point. Once p develops, it drives both S and 7\ 
order via the interactions 

f pQ = -Z WpQ plS (5.3) 
f P T = ^^JlwpQTplTi , (5.4) 

so that in mean-field theory [p4f 

S~ \r p ~r pc \ 2 ^\ (5.5a) 
Tx ~ |r p -r pc | 3 ^(3), (5.5b) 

for r p < r pc in the V phase. 
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2. I — > N transition 



,0(7) 
Jt 



u T \T\ 



(5.6) 



3D: 






top: 



FIG. 9. A schematic for the I — > N transition. 

The I — N transition is driven by the development of 
Q % i order from the Isotropic phase. There are no cou- 
plings that explicitly drive either p l or T^ k order once 
Q y order develops. Consequently, this well-studied tran- 
sition pl| , ^2| . which is described completely by the /q 
part of the free energy density, is generically first-order 
and in mean-field theory takes place at tq = Wq/12uq. 
A direct transition from the isotropic to the biaxial ne- 
matic (N + 2) phase is also possible. Since it is fairly 
complicated and has been treated in detail we will 
not consider it further here. 



where T is a 7 -dime nsional vector with components T M 
defined by Eq. ( 2.17c ). Because the underlying 0(3) sym- 
metry of our system is lower than 0(7) , with T ljfc forming 
its 7-dimcnsional irreducible representation, it is not sur- 
prising that the full free energy is not 0(7) invariant. The 
vt quartic term, explicitly breaks the 0(7) symmetry and 
determines which of the 7 T M irreducible components of 
T Ijfe order at the transition from the Isotropic phase. 

To determine which components order, it is convenient 
t o use the alternative representation of T % i k given by Eq. 
( p. 19 ) . In this representation we have 



/t = ifTT 2 



u T T A + f v 



(5.7) 



with 

JVT 



[9 cos 4 0i + ^ sin 2 20i (1 + 2 cos 20 2 ) 
25 sin 4 0i sin 4 2 



(5.8) 



breaking the 0(7) invariance of fx- fv T (and thus /t) 
has an 0(2) invariance, whi ch, th rough our choice of pa- 
rameterization of T yfe , Eq. ( 2.19| ), is manifested by /„ T 's 
being independent of 63. Thus, finding the minimum- 
energy state requires the minimization of /t(0i, 02 ) over 
two rather than three angles. 



3. Transitions from I driven by T % i 



(i.) I — > N + 3 transition 



Transitions from the / phase involving the develop- 
ment of third-rank tensor (2^j2^| order are more com- 
plex than the other transitions from the / phase we have 
considered. The many degrees of freedom in the T^ k 
tensor lead to the possibility of two distinct transitions, 
the / -> T and the I -> + 3 transition (27) . Since it 
is fluctuations in T yfc that drive these transitions, the 
noncritical degrees of freedom, p l and Q JJ can safely be 
integrated out to produce an effective theory involving 
only T 1 - 7 whose free energy is identical in form to Eq. 
(4.4). This is the energy that we will use to study tran- 
sitions from the I phase involving T^ k order. We will, 
however, have cause to return to the more general theory 
in our discussion of the / — > A^ + 3 transition. There are 
two im port ant things to note about the free energy fx 
in Eq. (14). First, in contrast to /q, this energy has no 
odd-order invariants because none can be formed with a 
third-rank tensor, T 1 - 7 . Second, there are two fourth or- 
der invariants ]27| ], which as we shall see, compete in the 
determination of the symmetry of the order parameter 
that develops from the Isotropic phase. In the limit of 
vanishing vt coupling, fx is invariant under the opera- 
tions of the group 0(7), as can be seen by reexpressing 
/t with vt = as 



3D: 




top: 





FIG. 10. A schematic for the / — » N + 3 transition. 

Following standard analysis, for vt > we find that 
the global minimum of /„ T (0i, 02 ) at fixed T is given by 



7T 

2 ' 
0, 



(5.9a) 
(5.9b) 



which corresponds to a state with 



11 



T 



T 2 IT 



(5.10) 



i.e., a s tate with a planar triadic order [see Eqs. (3.1) 
to (3.4)1, nere chosen to lie in the m — 1 plane. A ro- 
tation within this plane shows that for a more general 
choi ce [t han that defined by the representation choice, 
Eq.(pl9|)] ofmandl axes relative to the molecular body 
axes, such triadic planar order is described by an arbi- 



When vt < the global minimum of f VT (&i,02) at 
fixed T is given by 



^min 



7T 

2 

7T 



trary linear combination of the and T^ K tensors, cor- which corresponds to the state with 
responding to nonvanishing T2 and T3 order parameters, 
i.e., a nonvanishing T^. 

Since this triadic order defines a plane that brings with 
it a normal invariant under reflection, it necessarily in- 
duces uniaxial nematic order, Q y = 5(n*n? — \fi %3 ) with 
n along the normal. To lowest order, the develo pmen t of 
S is brought about by the Jqt coupling of Eq. ( 4.6c ), 



T ijk 



rpijk _ rri rijk 



(5.15a) 
(5.15b) 

(5.16) 



1 



/qt = ^w QT T 2 S , 



(5.11) 



which leads to the expected uniaxial nematic order with 



S = 



T 2 
fg 

2r Q 



(5.12) 



Hence the state for vt > is the N + 3 phase in which 
the nematic and triadic order, transverse to the nematic 
axis, coexist. From the point of view of symmetry it is 
equivalent to a liquid of orientationally ordered equilat- 
eral triangles with aligned normals. 

Another solution that minimizes the energy /t(^i, O2) 
and that is deg ener ate with the state described by the 
solution in Eq. (5.9) is 



^min 



arccos - 
n 

2 ' 



It corresponds to a state 
T 

71 



(5.13a) 
(5.13b) 



(5.14) 



which is equivalent to T^ ia = TI^ k , Eq. ( ^3 .10 ) after n 



and m are interchanged. Clearly then this solution also 
represents the N + 3 phase, but with the nematic axis 
along m rather than n, and the triadic order in the n — 1 
plane. 

(ii.) I —* T transition 



which is invariant under the operations of the tetrahe- 
dral group T<j. As can be seen in Fig. ^, the group Td has 
three Ci axes coinciding with the axes of the cube, four 
C3 axes coinciding with the body diagonals of the cube, 
six reflection planes passing through each edge and bi- 
secting the opposite edge of the tetrahedron, and four S4 
improper rotation axes corresponding to the axes bisect- 
ing (four) sets of two opposite edges of the tetrahedron. 
Because this state lacks n — > — n symmetry, no nematic 
or any o ther order is induced by the coupling free energy 
/ c , Eq. ( |4.6| ). Because of its tetrahedral symmetry and 
because only the Tq component of T l ^ k is nonzero, we 
identify this state with the T phase illustrated schemat- 
ically in Fig. ||(a). Since f VT is independent of 6*3, an 
arbitrary linear combination of Tq and T7, rather than 
Tq alone, will in general become nonzero at the I — ► T 
transition. 




Tetrahedratic 
(T) 




3D: 





FIG. 11. A schematic for the / — > T transition. 



FIG. 12. A portion of a phase diagram for banana-shaped 
liquid crystal, illustrating two possible transitions out of the 
Isotropic phase. For vt > 0, the transition is I — > iV + 3, and 
for vt < it is the I — > T transition. Although in mean-field 
theory these transitions are continuous, we expect thermal 
fluctuations to drive them first-order. Lowering temperature 
along a finely-tuned vt — curve, we expect a continuous 
transition in the 0(7) universality class. 
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For vt < we also find another solution that mini- 
mize s the energy and is degenerate with the state in 
Eq. (5.15). It is given by 



can be easily evaluated. The harmonic free energy / W 
can be expressed as a sum of five independent parts: 



(AO 



^min 
^min 




and corresponds to a state with 



(5.17a) 
(5.17b) 



(5.18) 



/ 



where 



?(AT) 



(AT) 



fW _ lf(AT)„2 , l f (AT) T 2 , ~W T C^oij,^ 



l^(Af) 



However, it can easily be shown that this solution is 
equivalent to T^ n 



2" r T 2|3 



Tit > after a rotation around the 1 Tl 5 = |fW (p3 + p |) + | f g0 ( T | + T 5 2 ) 



axis: 




(5.19a) 
(5.19b) 



It, thus, also corresponds to the T-phase with pure T§ 
order in a rotated coordinate system. The corresponding 
phase diagram that graphically summarizes phase tran- 
sitions outlined above is given in Fig. Il3. 

For vt = 0, it is clear from Eqs. ( |5.7| ) and ( |5.8D that 
the transition from the / phase is in the 0(7) universal- 
ity class, and it is to a state that spontaneously breaks 
0(7) symmetry by picking out a particular direction in 
the 0(7) symmetric space (a point on a 7-dimensional 
sphere) for the vector Ij, to point in. Clearly, as we have 
seen above, the vt coupling is relevant in the ordered 
phase and drives the resulting state toward N + 3 phase 
for vt > and toward T for vt < 0. 

We have also investigated the stability of the 0(7) sym- 
metric transition (with vt — 0) to a finite value of the 
vt symmetry breaking interaction |2{|. In a renormal- 
ization group calculation, just below the upper-critical 
dimension d = 4, we find that in the presence of thermal 
fluctuations the vt coupling always drives this transition 
first-order. This is analogous to the similar phenomena 
known in magnetic systems in which cubic crystal fields 
drive the 0(3) transition of hypothetical isotropic mag- 
nets first-order 32 33] . 



B. Transitions from the TV phase 

As illustrated in the flowchart of Fig. [?|, there are 5 
symmetry-reducing transitions from the nematic phase. 
These are the N—>V,N—>N + 2,N—>N + 3, 
N — > V + 2, and N — > Nt transitions, all of which we will 
discuss in detail below. To determine which transitions 
will occur for a given set of phenomenological parame- 
ters, we focus on the part of the full free energy density, 
f( N \ describing harmonic fluctuations about the nematic 
phase with nonvanishing S. This free energy, f^ N \ is de- 
termined by the harmonic parts of the energies f p , /q, 
and /t, and by the coupling terms f c . The most impor- 
tant contributions to f c come from f p Q, and / p qt, which 



a 



Pi,2,T 4 , 5 (Pl T 4+P2T 5 ) 



/(AO _ ld N )/n-2 , t2 



: -^;>{Ti+Tt) 



(5.21b) 
(5.21c) 

(5.21d) 
(5.21e) 



with 



r (AO = 
pa 

~(A0 



*P 3 w pQ^, 



r T - j^wqtS - (j^wi + 5W2 + gW 3 )S 



~(A0 

a P3,Ti 



sWpQtS, 



^ 2 =2r Q + Aw Q S+fu Q S', 



T T 2 ,3 = r T+ ^WqtS - g(u>i + W 3 )S , 



~(N) _ 2 
' Pi,2 — 1 P ^ 3 



w p qS, 



• 45 



Wi 



15 



W 2 



45 



w 3 )S 2 



iv 2 ry 



; = vt ~ f(2 Wl - w 3 )S z 



(5.22) 



Within mean-field theory, the N phase becomes unstable 
to the development of biaxial order (characterized by a 
linear combination of B\ and B2), of triaxial order (char- 
acterized by a linear combination of Ti and T3), and of 
Nt order (characterized by a linear combination of Tq 



and T7), when fg 3 |36|, ? T J^, and f Tfj7 , respectively, 
pass through zero. The N phase becomes unstable to lon- 
gitudinal vector order (characterized by coupled p 3 — T\ 
order parameters) and to transverse vector order (char- 
acterized by coupled p\ — T4 order parameters) , when the 
determinants, 



AN) 
X P3,Ti 



,(A0 



'Ti 



r Pl,2 r T 4 , 5 



(a [N) V 



(5.23a) 
(5.23b) 



respectively, pass through zero. The nature of transitions 
out of the N phase will be determined by which mem- 
beroftheset = I A (iv) ? {N) ? {N) f (N) A {N) i 
first passes through zero on lowering the temperature. 
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1. N — ► V transition 



2. N -> iV + 3 transition 



3D: 



top: 









FIG. 13. A schematic for the N —* V transition. 



The TV — > V transition is signaled by the development 
of vector order along the unique direction n of the N 
phase, i.e., by the development of p^ and T\ order. Thus, 
this transition occurs if A„ ^ is the first of the set Sm 



to become zero. The relative sign of p% and T\ is fixed 
by the eigenfunction associated with the smallest eigen- 
value of the matrix defined by fp^j< 1 ■ The overall sign is, 
however, arbitrary. Since it is the associated Zi symme- 
try (with V pointing along or antiparallel to n) that is 
broken, the N — ► V transition is in the well-studied Ising 
universality class with coupled order parameters 



P3 



(JV) 

P3,Tl 

(N) 

P3,Tl 



(5.24a) 
(5.24b) 



growing for A 



(N) 
P3,Ti 



< A 



(iV)c 

p 3 ,Tr 



Once and T\ become 



nonzero, they do not force the development of any other 
order, and the V phase is completely characterized by p^, 
S, and T\ order parameters, as illustrated in Figs. [| and 



3D: 



top: 





FIG. 14. A schematic for the N -> N + 3 transition. 

The N + 3 phase develops out of the N phase with 
the appearance of a linear com bin atio n of Ti and T3 or- 
der. As discussed in Eqs. (fTll) to (|3~^ ), T 2 and T 3 define 
a two-dimensional representation of the group of rota- 
tions perpendicular to n and describe triadic order in 
the plane perpendicular to n. Since the N phase is in- 
variant with respect to arbitrary rotations about n, the 
free energy of the N phase is a function only of the rota- 
tionally invariant combinations Tf + T$ = | , 3 1 2 - Thus, 
the N — > N + 3 transition is in the well-known XY uni- 
versality class. Within mean-field theory, this transition 
occurs when 7j? 3 is the first in the set Sn to pass through 
zero on cooling. T2.3 order drives no other order, and the 

N + 3 phase is completely characterized by S and T2,3 
with 



T 2 oc I3 



\f {N) 



r T 2 , 3 c\ 



(5.25) 



for fTn 3 < f. 



'T 2i3 c- 



3. N -> N + 1 (= V + 2) transition 



3D: 



top: 
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FIG. 15. A schematic for the N — » N + 1 transition. 



4- N — > N + 2 transition 



The V + 2 phase is distinguished from the N phase 
by the existence of a vector order described by p in the 
plane perpendicular to n. We can arbitrarily choose p 
to be along m so that p\ is nonzero. Since there is al- 
ready uniaxial order in the N phase, T4 order has the 
same symmetry in the N phase as does p\ order, and 
not surprisingly p\ and T4 (as are generically pi and T5) 
are coupled in fW, Eq. ( gjgjg ). The invariance of the 
N phase with respect to arbitrary rotations about n im- 
plies that the energy of the N phase must be a function 
of rotationally invariant combinations p\ + p\ = \pi,2\ 2 , 



T 2 



n 



|T 4j5 | 2 and piT 4 



P2T5 = pi, 2 • T 4i 5 as the 
harmonic free energy j^f^Tks in Eq. ( 5.21d| ) is. The 
N — > V + 2 transition is thus in the XY universality 



(N) 

Pl,2,T 4 , B 



IS 



class. It occurs in mean-field theory when A 
the first of the set Sn to pass through zero upon cooling. 
The order parameters p\ and T4 (and pi and T5 related 
to them by a rotation in the m — I plane) will drive the 
N — > V + 2 transition, both growing continuously from 
zero as 



Pi oc T 4 



v (iV) 



(5.26) 



for A^y 4 < A^"^. Once these order parameters be- 



(N)c 



come nonzero, however, they pick out a direction in the 
plane perpendicular to n that drives the development of 
a non-vanishing biaxial order B\ via the f v n and Jqt 
coupling free energies, as is clear from Fig. |l5j. In mean- 
field theory, B\ ~ p\ and T2 ~ v\- Below the critical 
dimension d c — 4, however, potentials in the coupling 
energies are relevant, and p4|,p5| 



3D: 



top: 





FIG. 16. A schematic for the N — > N + 2 transition. 

The development of nonvanishing biaxial order param- 
eter i?i,2 = (BijBz) converts the uniaxial N phase to 
the N + 2 phase. The biaxial order parameter is a rank-2 
symmetric-traceless tensor, which, because of its confine- 
ment to the two-dimensional plane perpendicular to the 
uniaxial axis n, is equivalent to a complex order param- 
eter forming an irreducible L = 2 representation of the 
U(l) group. The N — > N + 2 transition is thus in the 
XY universality class. In mean-field theory, it takes place 
when rs-12 * s the first of the set Sn to pass through zero 
and more generically, in the presence of thermal fluctua- 
tions, we expect, 



B x cx B 2 
for fgj < f W 



l r Bi,2 



'Bl,2cl 



(5.28) 



Biaxial order forces no p l or T y '* or- 
der, so the N + 2 phase is fully characterized by nematic 
order parameter S and an arbitrary linear combination 
of B\ and Bi biaxial order parameters, i.e., by -Bi,2- 



To 



Pl ! 



(5.27) 



5. N — > Nt transition 



where a n - 
dent on n. 



■ x n n[n — 1) with x n only weakly depen- 



To emphasize the secondary role that B\ and T2 or- 
der parameters play at the N — ► V + 2 transition, in 
Fig. (0), B\ and T2 are placed in parenthesis along the 
N — > V + 2 line. Therefore, the V + 2 phase is reached 
from the N phase (characterized by a finite value of S) 
via the well-studied, second order XY transition upon 
the development of p\ , B\ , T2 and T4 order parameters 
or a set that can be obtained from this one by an XY- 
rotation. 





FIG. 17. A schematic for the N — > Nt transition. 



15 



Finally, the Nt phase is distinguished from the N 
phase by the development of an arbitrary linear combi- 
nation of the Tq and T7 i.e., of the Tq^ order. Since such 
an order parameter picks out a single direction within 
the isotropic plane perpendicular to n, this transition, 
like the N -> N + 2, TV -> N + 3, and N -> V + 2 
transitions, is in the XY universality class. In order for 
this transition to occur, in mean-field theory rT e 7 must 
pass through zero before any of the other members of 
the set iS/v- If we restrict the interaction energies to / p q, 

for and / p qt, will never be the smallest in the set. 
However, higher-order terms of the form Jq^t 2 can favor 
the formation of T e — T7 order over the others and make 
this transition possible, with 

T 6 ocT 7 ~|fW-fS/-, (5.29) 

for < „. Because the of Tr — T7 order drives no 
other order, the Nt phase is characterized by nonvanish- 
ing S, and T 6 7 order. 

C. Transitions from the V phase 



'12 =42 

AV) _ ~{N) „ R 2 

a pi,2,T 4 , 5 - a pi,2,T 4 , 5 " Z \IZ W PTP^ 

U i n ) rr = ^WQTTi ~ 

Q3,4,Xi,e 5^3 w 1 V 15 

-^73(^1+ 6^2+4^3)5^, (5.32) 

As Fig. |?] indicates, there are symmetry-lowering transi- 
tions from the V phase to the V + 2, V + 3, and iV + V 
phases. Which one occurs is determined by which the set 

Sy = {f£>£U.r><U.W> where 

A« T =f { P P T V) - T )\ (5.33a) 

(and q 3 4 j- 4 r unenlightenly complicated) first 

reaches its critical value. The nature of the transition 
out of the V phase will be determined by which of these 
three mean-field reduced temperatures first reaches the 
respective critical temperature. 



In the V phase, three order parameters, ps, S, and T\ 
are nonzero. Harmonic fluctuations of the other order 
parameters in this phase are described by the free energy 
density 

f<y) = ?P T +/f)+/f) T , (5.30) 
where 

/£,U, 7 = ¥ { Z^l + & + |C(T| + T-) 
+&%] T6 (B 1 T 6 + B 2 T 7 ) 

C = lC(^ + r 3 2 ) (5.31) 

S,Q M ,r 4 , 5 = WM +*>!) + ¥Z]M + Ql) 

( T 4 + if) + 52,g, i4 (pi03 + P2Q4) 

+«pIU>i t 4 +p 2 T 5 ) + 5^ ) 4iT4i5 (Q 3 r4 + Q 4 T 5 ) 

where 

<l=<l -§(2»i +«*)!?, 

Q Bi,2,T 6 , 7 = TU^O^l " \J3 W PQTP3- 



1. V — > F + 3 transition 




FIG. 18. A schematic for the V — » V + 3 transition. 
This transition occurs if reaches its critical value 

^2,3 

f £c before A b^2,t 6 , 7 and A J,l2,Q3,4,T 4 ,5 reach their re- 
spective critical values, both zero in mean-field theory. 
Since there is rotational invariance in the space defined 
by T 2 and T 3 , this transition is in the XY universality 
class and corresponds to development of triaxial order 
T 2i 3 in the plane perpendicular to the vector order axis 
P3. Since no other order is driven by this development 
of the T 2)3 order, this order parameter, together with p 3 , 
S and T\ (already present in the V phase) completely 
characterizes the V + 3 phase. 



1G 



2. V — ► V + 2 transition 



3D: 





Thus the new phase has C\h symmetry and is the N + V 
phase. Once these order parameters have been estab- 
lished, they will drive B\, T2, and Tq nonzero via higher- 
order terms such as f p Q and f p r- Since the V — > N + V 
transition is controlled by the development of vector or- 
der in a plane, we expect it to be in the XY universality 
class. 



D. Transitions from the N + 2 phase 



top: 




FIG. 19. A schematic for the V — > V + 2 transition. 



If aw t 

01,2,16,' 



is the first in the set Sy to reach its crit- 
ical, there will be a transition from the V to the V + 2 
phase signaled by the development of biaxial order in the 
plane perpendicular to n characterized by a specific lin- 
ear combination of B12 and Tq^. Rotational invariance 
in this plane implies that the transition will also be in 
the XY universality class in which an arbitrary linear 
combination of (B\,Tq) and (B2, Tj) will spontaneously 
order. 



3. V — > N + V transition 



3D: 





top: 




FIG. 20. A schematic for the V -> N + V transition. 



If A 



(V) 



is the first in the set Sv to reach its 



A rp IB ill Lilt Sdl JV 

Pi, 2, V3, 4, -14,5 * 

critical value, there will be a transition from the V phase 
to one in which vector (^1,2) and third-rank tensor (T45) 
order develop in the plane perpendicular to the already 
existing vector order, p 3 . Without loss of generality we 
represent this transverse vector ordering by nonzero p\ 
and T4. This set of order parameters is invariant under 
only one non-trivial operation: reflection from y to —y. 



The N + 2 phase is the standard biaxial nematic phase 
with nonvanishing S and an arbitrary linear combination 
of biaxial B\ and B2 order parameters; without loss of 
generality, for convenience, we choose our coordinate sys- 
tem (m, 1, n) so that the biaxial order is described by B\. 
As the flow chart of Fig. Q indicates, the N + 2 phase can 
undergo symmetry-lowering transitions to the V + 2 and 
(Nt + 2)* phases. Harmonic fluctuations about the iV + 2 
phase, which determine the nature of transitions out of 
the N + 2 phase, are described by 



/ 



(N+2) _ f(N+2) 
- J P 3,T U T 6 



f(N+2) 
Jpi,T2,T 4 



f(N+2) f(N+2) 
J P2-T 3 ,T 5 ' JT 7 



(5.34) 



where 



;(JV+2) 

P3,T U T 6 



l~(N+2) n 2 , 1~{N+2) T 2 , 1 ~(N+2) T 2 
2 ' P3 ^3 + 2 ' T, 1 1 2 ' T 6 J 6 



. -(JV+2)^, ^ 



2'Ti 

• "p^rf^sTi + a p ^ 6 2) p 3 T 6 , (5.35a) 



?(JV+2) 
J Pi,T 2 ,T 4 



2>i Pi 



l-(iV+2) 

2'T; 



1 -(^+2)^2 



2 r T, 



/- 



i T2 .T i 'T2T A + & p N i 1 r l ) piT2 

'(N+2) _ l~(N+2) r 



- a 



] piT 4 



)rp2 



(5.35b) 
(5.35c) 



where 

? (JV+2 



'PS 



P3 

-(iV+2) _ ~(JV) 

'Ti ~~ 'Ti 

: (JV+2) _ _(JV) 

T 6 _ r T 6 , 7 

(iV+2) 
: Ti,T 6 



O: 



-1(2^+^3)^, 

-l{2 Wl + 2w2+w 3 )Bj, 
= ^k w QTBi-^ m (2w 1 +3w2 



a 



(N+2) _ _(iV) 
'P3,Ti -" 



! P3,T!' 



a 



(N+2) 
■P3,T 6 



j:(N+2) = ~(N) 



'in 



Pi, 2 



T*2 ^2 , 3 

_(iV+2) = ~(JV) 

T4 T4 , 5 



-W+2) _ 2 



X T 2 ,Ti 



(N+2) 



a 



Pi,T: 



■WpQrBi, 
2w pQ Bi, 
2( Wl +w 2 )Bf, 
lw QT B 1 + ^(w 1 (iS-7B 1 ) 
w 2 {AS - B x ) - 2w 3 {2S + Bx))B x , 

+3w 2 (2S - B x ) + w 3 (2S - ZB X )B X ), 
= —WpQxBi, 
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-(JV+2) _ ~(JV) 



-(JV+2) _ ~(JV) 
r T 7 - r T eJ 



\(2w 1 -w 3 )Bl 



(5.36) 



and fp^T^Tr IS obtained from J^t^Ta D ^ replacements: 
Pi -> P2, r 2 -> T 3 , T 4 -> T 5) Bi -> -Bi. 

Following our earlier analysis of other transitions, 



we introduce the determinants A 

(N+2 



(N+2) 



(N+2) 



,Ti,T B > pi,T 2 ,T 4 



and 



P2 t 3 Tr °f the matrices, respectively defined by the co- 
efficients in the free energy densities fp^Tx,T e i fpuTab 



and f {N+2) 



Which transition from the N + 
phase takes place is determined by the set iSat+2 



{r 



~(N+2) A (N+2) 



,A 



(N+2) 



P3 ,T X ,T 6 J " Pl ,T 2 ,T 4 ' ^p 2 ,T 3 ,r t 



A 



(JV+2) 



}• 



i. A?" + 2 ^ F + 2 transition 



3D: 






top: 



FIG. 21. A schematic for the N + 2 — > V + 2 transition. 

Interestingly, there are three routes from the iV + 2 
phase to a phase with V + 2 symmetry. In the first route, 
the nematic axis along n, which in the N + 2 phase is 
invariant under n — > — n, is converted to a vector axis 
with the development of p 3l T±, and Tg order. The bi- 
axial order of the N + 2 phase persists, resulting in the 
V + 2 phase with vector order along the two-fold axis n. 
In mean-field theory, this transition to V + 2 phase takes 

place when A p ^^T 6 ^ s the tirst m s °t &N+2 to pass 
through zero. 

In the second route, vector order develops along the di- 
rection m, corresponding to the maximum eigenvalue of 
the nematic order parameter Q lJ in the plane perpendicu- 
lar to n. In this route, which, in mean-field theory occurs 



when A 



(N+2) 
Pi, t 2 ,t 4 



is the first in Sn+2 to go through zero, 



the pi, Ta, and T4 order parameters become nonzero. 

In the third route, A^^ 2 ^ is the first in Sjv+2 to go 
through zero, and vector order develops along the 1 axis 
perpendicular to m, defined by the eigenvector of Q lJ 
with minimum eigenvalue in the plane perpendicular to 
n. It is the P2, I3, and T5 order parameters that become 
nonzero at the transition. 



In all three versions of the N + 2 — > V + 2 transition 
the vector p order develops along one of the three two- 
fold symmetry axes defined by the biaxial order of the 
N + 2 phase, the eigenvectors of Q 1 - 7 . In each case, p can 
develop either a positive or negative value along an axis 
already chosen by the N + 2 phase. Consequently, these 
N + 2 — > V + 2 transitions are in the Ising universality 
class, with well-known critical properties. 



2. N + 2 -> (N T + 2)* transition 



3D: 





top: 



FIG. 22. A schematic for the N + 2 — > (N T + 2)* transition. 

For our choice of the m-l axes, for which the biaxial 
order of the N + 2 phase is characterized by Bi ^ and 
B2 = 0, another symmetry-reducing transition from the 
N + 2 phase is signaled by the appearance of the T7 or- 
der. It corresponds to the N + 2 — > (Nt + 2)* and takes 
place when the reduced temperature f^! V+2 ' ) is the small- 
est in the set iSjv+2. For a different choice of axes, in 
which B2 ^ and B\ = characterize biaxial order, this 
transition would instead correspond to development of 
Tq order. And, more generally, in an N + 2 phase char- 
acterized by a particular linear combination of B\ and 
B2, i.e., by B12 a corresponding orthogonal linear com- 
bination of Tq and T7, i.e., the T§j order parameter, such 

that Bi i2 ■ f 6 ,7 = develops at the N + 2 -> (Nt + 2)* 
transition. 

To understand the nature of this transition it is helpful 
to first visualize the N + 2 phase as a collection of ori- 
cntationally ordered but positionally disordered planar 
diamond units, each consisting of two, leg-to-leg banana 
molecules, illustrated in Figs. §(b), ||(c), and [22]. One 
can choose the {n, m, 1} triad such that the diamonds 
lie in the n-m plane and, therefore, that their biaxiality 
is characterized by nonvanishing S and B\ order param- 
eters. Now the transition to the (Nt + 2)* phase cor- 
responds to counter, out-of-plane twisting about the n 
axis of the two diamond-forming banana molecules. It is 
signaled by the development of the T7 order parameter, 
with the twist angle 5 in the range < 6 < tt/2 given by 
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tan (5 = T7/B1. 



(5.37) 



We note that the lower limit of this range 6 = cor- 
responds to the planar diamonds of the achiral N + 2 
phase and the upper limit, <5 = tt/2 (reached only in the 
limit Tj — ► 00 or B\ = 0) corresponds to the Nt phase, 
which is also achiral, and lacks biaxial order . In con- 
trast all other values of the twist angle < 5 < tt/2, 
corresponding to a nonzero value of T7 describe a spon- 
taneously induced chirality of the (Nt + 2)* phase. Since 
the sign of the twist angle (and correspondingly the sign 
of T-j) can be either positive (right-handed) or negative 
(left-handed), in mean-field theory the N+2 -> (N T +2)* 
transition is clearly in the well-known Ising universality 
class (but see below). 

Finally, we observe that in a phase which sponta- 
neously breaks chiral (mirror) symmetry, a tensor, totally 
antisymmetric in all its indices (akin to the well-known 
tijk tensor) must spontaneously emerge. It is easy to 
verify that in the chiral (Nt + 2)* phase, most gener- 
ally characterized by finite S, B\, B2, Tq and T7 order 
parameters, is given by 

E ijk = Q™B 3t T stk + Qi s B kt T sti + Q ks B lt T stj , (5.38a) 
= j= SfaT-r - B 2 T 6 ) ei jk , (5.38b) 

is indeed such a fully antisymmetric tensor, which devel- 
ops spontaneously from our theory, that is based solely 
on fully symmetric tensors Q y and T 1 ^ . 

The existence of a fully anti-symmetric third-rank ten- 
sor allows for invariant couplings linear in spatial deriva- 
tive. For example, a term of the form 



E ijk Q il VjQ kl ~ T 7 n • V x n 



(5.39) 



is permitted, where the right hand side represents the 
dominant part near the transition where S and B\ are 
effectively constant, T7 is small, and n is the Frank di- 
rector. The elastic energy of the N + 2 phase includes 
the usual twist energy K2 [n • V x n] 2 . There are other 
terms in the elastic free energy arising from the biaxial 
order of the N + 2 phase, and they may influence the 
nature of the ground state of the (N t + 2 )* phase. The 
existence of terms such as that of Eq. 5.39| that are linear 
in spatial gradients implies that the ground state of the 
(Nt + 2)* phase will be spatially inhomogeneous. The 
simplest chiral inhomogeneous phase we can imagine is a 
cholesteric (Nt + 2)* phase in which n rotates in a helical 
fashion as in the standard cholesteric with pitch P and 
pitch wave number qo = 2tt/P which near the transition 
to the N + 2 phase at temperature T c scales as 



go 



T 7 /K 2 ~ I ATI 1 / 2 



(5.40) 



in mean- field theory where AT = T — T c . In a true 
critical theory, qo will also scale to zero as a power of 
I AT I provided the transition remains second order. In 
the presence of fluctuations, the T711 • V x n term is likely 



to modify the universality class of this transition to some- 
thing other than the naively expected Ising universality 
class (possibly even driving the transition first-order), 
but we have not, however, analyzed the critical theory 
in detail. If molecules themselves are chiral, or a chiral 
dopant is added, then the N + 2 phase will be a chiral 
(N + 2)* phase with a nonvanishing qo. Chirality will act 
like an external ordering field (like the magnetic field of 
an Ising model) for qo with strength h, and if the N + 2- 
to- (Nt + 2)* transition is second order, one can expect 
qo (T, h) to scale as 



q (T, h) = \T-T c f*f(h/\T-T c \ A ), 



(5.41) 



where T c is the transition temperature, A is the external- 
field gap exponent, and /3 q is a critical exponent, which 
in mean-field theory is the order-parameter exponent 
(3 = 1/2. 



E. Transitions from the N + 3 phase 

As we indicate in Fig. |?], the N + 3 phase, described 
by coexistence of the uniaxial order parameter S and the 
transverse (to the nematic axis) triaxial order parameter 
T2.3 = (T2,T^), can undergo symmetry- lowering transi- 
tions to the V + 2, V + 3, and N + V phases. 

These transitions are all signaled by the development 
of vector order p. They are distinguished by whether this 
vector order is along (N + 3 — > V + 3 transition), trans- 
verse to (N + 3 ^ V" + 2 (= N + 1) transition), or at an 
arbitrary polar angle < 5 < tt/2 to the uniaxial n axis, 
singled out in the nematic phase. Although it is conve- 
nient to think of these transitions as driven by the vector 
order p, more precisely they are driven by specific linear 
combinations of the vector and other order parameter, 
linearly coupled to p. 

It is convenient to choose the m-1 axes so that the tri- 
axial order of the N + 3 phase, in the plane perpendicular 
to n phase is described by T 2 ^ and T3 = 0. For this 
choice, the harmonic fluctuations about the N + 3 phase 
are described by 



f(N+3) _ f(N+3) f(N+3) f(N+3) 
■> ~ J P3,Ti ^ J Pl ,B u T 4 "T" JQ 3 ,4,T 6|7 i 



(5.42) 



where 



f(N+3) _ 1 ~(JV+3) 2 , 1 ~(N+3) T 2 , 5^+3)- 
J P3,Ti — 2 r P3 P3 + 2 r T ± J l T a p 3 ,Ti P3J-1 



fiX,T,. - m 3) & +p*) + ^(bi + BD 

+C!S>i T 4 +P2T5) + 5^,(5^4 - B 2 T 5 ), 

Cr - ^£: 3) («3 + Ql) + *CCtf + T?) 

+«£,t 3 T 6 , 7 (Q3T 6 -Q4T 7 ), (5.43) 
with the coefficients 
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f(JV+3) 

P3 

-(JV+3) 



~(JV+3) _ -(JV) 



f(JV) 

P3 ' 

^ N ^(Au T + lv T )Tl 



2. N + 3 -^V + 2 f= AT + 1 j transition 



a P3,Ti 
f(JV+3) 

Pl,2 

f (JV+3) 
-(JV+3) 

^4,5 

^(JV+3) 
(JV+3) 

Pl,2,T 4 ,5 

-(JV+3) 

£>1,2 ,-^4,5 

-(JV+3) 

r <53,4 



a 



f(JV) 

Pi, 2 ' 



rW-2( Wl+W2 )T 2 2 , 
fW+(4 UT + | UT )T|, 



-w pQT T 2 , 



Pl,2:T 4 , 5 

r Q - w Q 5 + |w Q 5 2 - iwiTf , 



^(^+3) _ -(JV) 
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75WqtT 2 - ^(6^1 - w 3 )^T 2 . (5.44) 



Transitions out of the N + 3 phase are controlled by the 
set«S JV+3 = {Ai JV + 33 ! Af + | rp ,A[f +3 » } of deter- 

iv^o L P3.-I1 ' Pi, 2, -Dl. 2. -14, 5 ' V3,4,lS,7 J 

minants of the harmonic coefficients that can be read off 
from /( Ar + 3 ) above. 



1. JV + 3— +V + 3 transition 



3D: 



top: 






FIG. 23. A schematic for the N + 3 — > V + 3 transition. 

The jV + 3 — > + 3 transition occurs when A Pl ; t 3 is 
smallest in the set Sn+z- In this transition, vector order 
develops along the nematic n axis to produce a linear 
combination of p% and T% order, with the latter explicitly 
induced as a third harmonic of the p^ order. The discrete, 
Z%, n — > — n symmetry, characterizing the jV + 3 phase is 
lost at this transition. Consequently, the N + 3 — * V + 3 
transition is in the familiar Ising universality class. 



3D: 





top: 



FIG. 24. A schematic for the N + 3 — > N + 1 transition. 

V + 2 order develops from N + 3 order by sponta- 
neously favoring one of the three equivalent directions 
in the plane perpendicular to n and, thereby, converting 
transverse triaxial order of N + 3 into transverse vec- 
tor order of N + 1, equivalently described by V + 2. 
With p chosen to be along m , the transverse vector or- 
der is signaled by the development of finite p\ order pa- 
rameter. As can be seen from f„ + S 

J Pl,2,-fc>l,2 



Eq. (5.43) 



order parameters B\ and T4, associated with p\ and lin- 
early coupled to it, are also simultaneously induced at 
the N + 3 — > V + 2 transition. For a more generic choice 
of the orientation of m-1 within the plane perpendicu- 
lar to n, the V + 2 phase is described by a specific lin- 
ear combination of p\^ = {pi,P2), -81,2 = (^1,-82), and 
T 4j5 = (T4, T5), obtained from the basic set {p\,B\,T±} 
by a planar rotation about the nematic axis n. 

The N + 3 — > V + 2 transition is the exit of choice out 
of the N + 3 phase, when A„ +3 J ™, is the smallest in 

r ' Pl, 2,^*1,2,-14,5 

the S/v+3 set. Since the N + 3 phase is characterized by 
three equivalent directions (in the plane perpendicular to 
the nematic axis n) along which vector order p can de- 
velop, we expect this transition to be in the universality 
class of the 3-states Potts model, believed to be weakly 
first-order in three dimensions and continuous in two. 

These symmetry based expectations are born out by 
our detailed computations, which show that the rota- 
tional degeneracy in the m-1 plane, which is present in 
the harmonic free energy, fp^Bi 2 t 4 5 > Eq. ( |5.43 ) is lifted 
by energy contributions of the form 



J pi 



'(JV+3) 
,-Bi, 



!,T 4 



= —a p COS 30 — OtB cos I 



(5.45) 



where <f> is the angle between the developing transverse 
vector order parameters pi, 2 (as well as B\ 2 and T^) 
and the m axis defined by triaxial order of the N + 3 
phase, 
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WpTPio 



5VT5 



v T T? 



4,.-, 



(5.46a) 
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UB = WQ3 T 2B 12 T 2 , 



and ^/pf^Bj 2 t 4 r arises from the following nonlinear 
couplings 

f P T = -w pT P Vp k T^ k (5.47a) 

/ Q3T 2 = -U; Q 3 T 2 QhhQi 2 frQhh T hi2i3 T hj2j3 ( 5 . 47 b) 

f T 4 — vt T ili2i3 T iliii5 T i2i4,i6 T i3i5i6 (5 47c) 

The three degenerate minima of the free energy 
^fp^Bi 2 t 4 r i Eq. (5.45), for the vector order parameter 
to settle into, correspond precisely to the three equiv- 
alent states of the 3-state Potts model, supporting our 
expectation that the N + 3 — > V + 2 transition is in the 
3-state Potts model universality class. 



(5.46b) 



a total of 6 free-energetically degenerate states. Integrat- 
ing out the Tqj order parameter and focusing on Q3 4 
alone, shows that the 6 states correspond to alignment of 
mQ along 3 vertices and 3 edge-bisectors of the equilat- 
eral triangle defined by the N + 3 phase. Thus there is 
a six-state clock symmetry described by a coupling pro- 
portional to T%Q% and we expect the N + 3 -> N + V 
transition to be in the universality class of the six-state 
clock model, which apart from irrelevant variables is 
in the universality class of the XY-model. Once Q3 
and Tq order is established, T4, pi, and B\ order is 
driven by couplings of the form Q\B\ and others. Thus, 
P\ ~ B\ ~ T4 ~ Q|. Then in mean- field theory, B\ or- 
der will drive P3 and Ti order via couplings of the form 
£iT 6 p 3 . Thus p 3 ~ T x ~ T 6 3 ~ Q|. 



3. A?" + 3 -> A^ + V transition 



F. Transitions from the T phase 



3D: 



top: 






The Tetrahedratic T phase with T4 symmetry is char- 
acterized by a nonvanishing arbitrary linear combination 
of Tq and T7, which we collectively refer to as Tq^. In 
this section, without lost of generality, we will choose the 
orientation of the m — 1 axes so that T7 is the only non- 
vanishing order parameter in the T phase. The harmonic 
free energy density for fluctuations from the T phase can 
be expressed as 



f(T) _ f(T) 7{T) fK , 
J ~ JQx + JQ2 + Jpi 



(T) 



(5.48) 



where 



FIG. 25. A schematic for the N + 3 — > N + V transition. 



The transition that occurs when A 



is the small- 



est of the set Sn+3 is, as we shall see, to the N + V phase, 
but the development of all of the order parameters char- 
acterizing this phase is complicated. At the harmonic 
level in f 6)7 and Q 3)4 , f (N+3) possesses a rotational in- 
variance in the plane perpendicular to the nematic axis 
n. However, because the N + 3 phase is characterized by 
3-fold order in the xy plane perpendicular to n (aligned 
along the z axis), this continuous symmetry will be re- 
duced to a discrete clock-symmetry. 

It is convenient to focus on ordering of Q3 and Tq. 
The established A^ + 3 order can be viewed as an equi- 
lateral triangle that favors alignment of m towards one 
of its vertices. Thus one might predict the symmetry 
of a 3-state clock model. However, the apolar (nematic 
along z axis) nature of the Af + 3 state ensures that the 
free energy f^ N+3 ^(Q3, Tq) is invariant under n — » n. 
This, combined with transformation properties of Q3 
and Tq under n — * n, with Q3,Tq — > — Q3, — T@ gau- 
rantees an additional Ising symmetry of the free energy 
/( w +3)(Q 3 ,r 6 ) = /(JV+3)(_Q 3; _T 6 ). This then leads to 



7(T) 
JQi 
f(X) 
JQ2 

f(T) 

'Pl,2,3<23,4,5 



if (T) 2 

2'Qi ^1' 
2 r Q 2 



Q21 



(5.49) 

rQ SA ,s (Ql + Ql + Qt) + \r p {p\ +p 2 2 + pi) 
^w pQT T 7 (p 1 Q 4 + P2Q3 + P3Q5) 



with the coefficients given by 



r Q 



r Q 3 



r Q - \{2wi 
r Q - \{2w\ 



™ 3 )T 7 2 , 

2w 2 + w 3 )Tj. 



(5.50a) 
(5.50b) 
(5.50c) 



The structure of /*- T ' implies that there are three 
symmetry-lowering transitions from the T phase, driven 
respectively by fluctuations in Q%, Q2, and a linear com- 
bination of pairs p\ and Q4, P2 and Q3, and P3 and Q--,. 
Which of these transitions occurs is determined by the 
set St -\f {T) r (T} A (T) i 

In the presence of T7 order, all of the p 1 ^ and Q?f or- 
der parameters exhibit third-order invariants. Since they 
play an important role in determining the nature of tran- 
sitions from the T phase, we display them here: 
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/q 3 



^mQ2(Q\ - Qi) - ~mQ\ ~ 75O3Q4Q5] 



2. T -> (JVt + 2)* transition 



fp*T = -V6w p tTtPiP2P3- 



(5.51) 



We note the appearance of the Q3Q4Q5 term in /q3, 
which can be paired with the P1P2P3 product in f p 3T- 
There are of course also third order terms in Q, whose 
coefficients are proportional to powers of T7, which, for 
small T7, are subdominant to wq term that we displayed 
above. Because these higher order terms do not qualita- 
tively change our results, with their effects accounted for 
by an effective wq coupling, we will not consider them 
here. 



1. T — » Nt transition 





FIG. 26. A schematic for the T — > Nt transition. 

As illustrated in Fig. 0, one of the three possible 
symmetry-reducing transitions out of the Tetrahedratic 
T phase is the T — * Nt phase transition. The Nt or- 
der develops by favoring one of the 3 two-fold axes of 
the tetrahedron through the growth of uniaxial Q % i or- 
der along that axis as shown in Figs. ||, |5[ and ^6|. In 
our parameterization, we focus on the two-fold axes de- 
fines by ±n. In this case, the T — ► Nt transition is 
signaled by the development of Q\ (or equivalently S) 
order. Becaus e of the existence of the third-order Q\ 
invariant, Eq. ( |5.5l| ), this transition is generically first- 
order. It occurs in mean-field theory when (T7) = 
7.2 



-w. 



(Tr)/12v,Q(Ty), where the third- and fourth-order 
couplings wq(Ty) and uq(Tj) can in principal depend on 
strength of the T order, characterized by T7 order pa- 
rameter. Because no other order parameters is explicitly 
induced at this transition, and with the convenient choice 
of the m-1 axes (that we have made here) the resulting 
Nt phase is fully characterized by the T7 and S order 
parameters. 





FIG. 27. A schematic for the T -> (Nt + 2)* transition. 

If, for a T phase characterized by the T7 order pa- 
rameter, Q2 (or equivalently B\) orders before Q\ or the 
Pi ~ Q3 combination, the transition from the T phase is to 
the spontaneously chiral (Nt + 2)* phase. If we focus on 
the m-1 plane, T7 order displays biaxiality (m l V + mH % ), 
with principle axes at 7r/4 relative to the biaxiality of the 
B\ order characterized by (vrtm? —1 % V). Since T7 and B\ 
are different order parameters, a simple rotation to de- 
fine a new biaxial axis is impossible, and the result is that 
the reflection invariance in the m-1 plane is spontaneously 
lost as B\ develops. The lack of reflection symmetry is 
the characteristic featur e of the chiral (Nt + 2)* phase, 
discussed in Sec. VD2. In analogy to the development 
of biaxial order from the Isotropic phase, once B\ orders, 
the Q1Q2 coupling in /q3 drives the development of S (or 
Qi) order. Thus, the (Nt + 2)* phase has nonvanishing 
S, B±, and T7 (or equivalently S, B2, and Tg) order. Al- 
though the T — > (Nt + 2)* transition is continuous within 
mean-field theory, based on the experience with the de- 
velopment of uni- and bi-axial orders from the Isotropic 
phase, we expect that here too, the T — * Nt transition 
will always preempt the T — * (Nt + 2)* transition. 

Within second-order mean- field theory, B\ ~ |AT| 1//2 , 
S - B\ - I AT I, and E ijk ~ T-jB^ - |AT| 3 / 2 . Since the 
T phase has elastic energies resisting spatial variations 
of the n, m, and 1 directions, the wavenumber of the 
cholesteric structure of the cholesteric (Nt + 2)* phase 
will scale as q$ ~ |AT| 3 / 2 in mean-field theory. 



3. T — ► V + 3 transition 





FIG. 28. A schematic for the T -> V + 3 transition. 
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The ordering of linear combinations of (pi,Qa), 
(P21Q3), an( i (P3,Q5), which we will refer to as p — Q 
ordering, leads to the V + 3 phase. This can be seen 
by ob serving that the third order potential fp*,T m Eq. 
(5.51) favors p\ — ±p 2 — ±J>3 with relative signs deter- 
mined by the sign of w p tT-j. Examination of Eq. ( 5.51 ) 
show that similar considerations apply to Q3, Q4 and 
Q5. Thus, the vector p will align along the (1,1,1) or 
a symmetry equivalent axes, i.e., along a three- fold axis 
of the tetrahedron characterizing the T phase. Thus, to 
discuss the phase transition signaled by the onset of p—Q 
order, its is useful to transform to a new coordinate sys- 
tem with n' along the (1, 1, 1) axis, which is most easily 
achieved by a rotation about the (1,-1,0) axis. Under 
this rotation we find that our basis transforms according 
to 



p(N T ) 



1(1 



1 



i = -3(i-7s) m ' + i(i 

n = * (m' + 1' + n'). 



73" 



73 1 



1 H' - 1 



V3' 



(5.52a) 
(5.52b) 

v aaa ~^ 1 1 xx j . (5.52c) 
Straightforward algebraic manipulations then yield 



T ijk _ V5 r/^fc 
h - T 1 1 



Tv J 2 + 1 3 )i 

so that in the rotated coordinate system, 



(5.53) 







(5.54a) 


n = 


-&T 7 , 


(5.54b) 


n= 


3 7 


(5.54c) 



In this rotated basis the polar order is described by 
p = p 3 n', also inducing the nematic order S' through 
the f p Q coupling. Thus, the phase produced by such 
p — Q ordering from the T phase has nonvanishing p 3 , 5", 
T[, T 2 , and T3 order, which describes polar order along 
the n' axis and triaxial order in the m-1 plane perpendic- 
ular to it; a rotation by 7r/12 about the (1, 1, 1) axis can 
be used to remove T3 component of the triaxial order. 
Thus, the phase produced by p — Q ordering inside the 
Tetrahedratic T phase is indeed the previously discussed 
V + 3 phase. 



G. Transitions from the Nt phase 

The Nt phase is characterized by nonvanishing ne- 
matic order parameter, S, and an arbitrary linear com- 
bination of the tetrahedratic order parameters, Tq and 
T7, which we collectively call Tq^. These define the di- 
rections of the orthonormal triad (m, 1, n) , which we can 
for convenience always pick to have T e vanish, with the 
tetrahedratic order completely characterized by the value 
of T 7 . 

The free energy f( NT ) describing harmonic fluctuations 
in the Nt phase can be expressed as a sum of three in- 
dependent parts: 



f(N T ) 
JBi 



f(N T ) , f(N T ) 
■ / P3,B 2 ,Ti ' J Pl ,2,Q 



(5.55) 



where 
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(5.57) 



As can be seen from the structure of f( NT > above, 
there are three possible symmetry-reducing transitions 
that can take place out of the Nt phase. In contrast to 
the nematic phase with Dooh symmetry in which fluctu- 
ations in the biaxial fields B\ and B 2 are degenerate, the 
Nt phase with nonvanishing Te,7 order breaks the de- 
generacy of fluctuations in B\ and B 2 ■ For our choice of 
m-1 axes, with Tq = 0, the tetrahedratic order parameter 
Tj couples the biaxial order parameter B 2 to the vector 
order along n, described by P3 and its third harmonic 
Ti. Hence two (of the three) transitions are the ordering 
of Bi, and the ordering of a linear combination of p 3 , 
B 2 and T\. A third possible transition out of the Nt 
phase is the development of vector order transverse to 
nematic axis n described by an arbitrary linear combi- 
nation of pi and p 2 . In the Nt phase with nonzero T7, 
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harmonic fluctuations in pi are coupled to those of Q4 
and T4, and harmonic fluctuations in P2 are coupled to 
those of Q3 and T 5 . As a result, the development of trans- 
verse vector order from the Nt phase is accompanied by 
a specific linear combination of these higher-order order 
parameters. 

Which of these three transitions takes place first is de- 
termined by minimum determinant in the set Sn t — 
{As, ,A^ T J rp , A„ T Vi rp rj, } of determinants of 

the harmonic coefficients, that can be read off from j^ T > 
above. 



1. Nt — > (Nt + 2)* transition 




FIG. 29. A schematic for the N T -» (A^t + 2)* transition. 

Development of biaxial order, Si, in the presence of T7 
(with Tq = 0) and 5 1 order is quite similar to the develop- 
ment of T-j order in the presence of B\ order, discussed in 
Sec. |VD 2\ and corresponds to spontaneous chiral sym- 
metry breaking of the Nt — > (Nt + 2)* transition. As can 
be seen from the structure of f^ NT \ no other order pa- 
rameters are induced at this transition and, as discussed 
above, the resulting (Nt + 2)* phase is characterized by 
nonzero S, B\ and T7 order parameters. Because it is the 
underlying Si — > —Si, Z2 symmetry of the Nt phase, 
that is broken when Si orders, the N T — ► (N T + 2)* 
transition is in the Ising universality class if the linear 
gradient coupling proportional to Sy^ is ignored. The 
latter term, aside from potentially modifying the critical 
properties of this transition, leads to a pitch wavenum- 
ber in the cholesteric (Nt + 2)* phase that scales in 
mean-field theory as go ~ B\ ~ IATI 1 / 2 just below the 
Nt -> (N T + 2)* transition. 



2. N T — > V + 2 transition 




FIG. 30. A schematic for the Nt — > V + 2 transition. 

The Nt — > V + 2 transition takes place via the devel- 
opment of vector order p = p 3 n along the nematically 
ordered n axis. Since in the presence of T7 and S such 
longitudinal vector order is coupled to S2 and Ti, the 
Nt — > V + 2 is also accompanied by the development of 
S2 and the T\ order parameter. Because the free energy 
of the Nt phase is invariant under n — » — n, the devel- 
opment of P3 longitudinal order can be of either sign and 
the transition is in the Ising universality class. 



3. N T -» N + V transition 

A third possible transition out of the Nt phase takes 
place when transverse vector order pi 2 develops. As can 

be seen from the form of fi^n T T , the develop- 

ment of p\ 2 is accompanied by biaxial order Q3.4 and 
a linear combination of T2 : 3 and T4 ; 5 order parameters. 
The resulting phase has aCn symmetry, and we there- 
fore identify it with the previously discussed N+V phase. 
At the harmonic level, the free energy appears to be 0(2) 
invariant, with respect to rotation of pi 2> However, in 
the presence of T7, nonlinearities in pi^, break this rota- 
tional invariance. A lowest order such symmetry break- 
ing nonlinearity is given by 

f N T _ ftiiM r rii,ji,kirpi 2 ,j2,k2 T iji r ik 1 j 2 k 2 

J symm. breaking — l * ± J V P V P i 

= ISTfrlpi (5.58) 

It introduces into the 0(2) symmetry harmonic free 
energy of pi ; 2 a well-studied cubic symmetry breaking 
anisotropy. Based on these studies [B2| , we therefore ex- 
pect the Nt — * N + V transition to be driven first-order 
by thermal fluctuations. fl33[ 
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H. Transitions from the V + 2 phase 

As can be seen from the flow-chart, Fig. fj], the biaxial 
vector V+2 phase, characterized for example by finite ps, 
S, B\, Ti, Tq order parameters and by Civ symmetry, can 
further lower its symmetry in two ways. It can undergo 
a transition to the N + V phase via the development of 
polar order, characterized by pi t 2, Qs,4, and T^, order 
parameters along one of the biaxial axis perpendicular to 
the V (j>s) order. Alternatively, it can undergo a tran- 
sition to the (Vt + 2)* phase via the development of T7 
order. Both the V+2 -> N+V and the V+2 -> (V T + 2)* 
are expected to be in the Ising universality class because 
in both cases it is Z2 symmetry that being broken. 



crystalline phases of bent-core molecules, particularly in 
phases with exotic symmetries not encountered in sys- 
tems of rod- or plate-like mesogens. Our work, however, 
suggests the possibility of smectic or layered analogs of 
these exotic fluid phases. In these putative phases, which 
we will explore in more detail in a separate publication 
Jl7t the smectic layer normal N provides an additional 
direction that may or may not coincide with a symmetry 
direction of the fluid phase. Most of the banana smectic 
phases classified to date are based upon V + 2 order em- 
bedded in smectic layers either in identical configurations 
(as in the SmCsPp phase in the notation j| of rcf. ||) 
neighboring layers or in alternating configurations (as in 
the SulCU-Pa phase) in neighboring layers. 



I. V + 3 — > N + V transition 

The three-fold symmetry in the plane transverse to the 
vector (p 3 ) axis of the V + 3 phase can be spontaneously 
broken with e.g., biaxial order in this plane driving the 
transition and other parameters (listed in Table 1) also 
condensing. We therefore expect the — > N+V tran- 
sition to be in the three-states Pott's model universality 
class. 



J. (JV T + 2)* -> {Vt + 2)* transition 




FIG. 31. A schematic for the (Nt + 2)* -» (Vt + 2)* tran- 
sition. 

The final transition that we will comment on is polar 
ordering transition from the nonpolar Nt + 2* phase. Be- 
cause of the present nematic order in Nt + 2* state, polar 
order breaks Z2 symmetry and we expect (Nt + 2)* — > 
(Vt + 2)* transition to be in the Ising universality class. 

VI. NEW SMECTIC PHASES 

Our primary interest in this paper has been in 
fluid (spatially homogeneous) but anisotropic liquid- 




FIG. 32. Schematic of smectic phases (a) SmAjv T , (b) 
SmAf, (N T + 2) , (c) SmAl (N T + 2) . 
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Each layer is then characterized by the molecular di- 
rections n, m, and 1 and by N. Though the V + 2 fluid 
phase is not chiral, each layer of smectic phases derived 
from it can be chiral if, for example, m lies in the smectic 
layer and n is tilted relative to N as is the case in the 
globally chiral SmC a Pa and SuiCsPf phases, and the 
SihC'sPa and SuiCaPf phases with chirality alternative 
in adjacent layers 

New types of smectics can arise from layering of Nt, 
(N T + 2)*, and (V T + 2)* fluid phases. In the simplest 
of these phases, which we label SuiAn t and depict one 
layer of in Fig. |32|(a), each layer has Nt symmetry with 
n parallel to N and m along a common direction in each 
layer. This phase has Did point group symmetry. Many 
variants of it are imaginable. For example, the m-axis 
could rotate by tt/2 from layer to layer, or the n axis 
could tilt relative to N to produce a SmC*(N T ) phase, 
which like the SmC* phase would be chiral. Layered 
phases formed from the (Nt + 2)* fluid phase are chiral. 
In one such phase, the Sirivlu (Nt + 2) phase, a layer of 
which is depicted in Fig. H^(b), the n axis is parallel to 
N, and the biaxial m axis rotates from layer to layer like 
the c director in the SmC* phase. If the pitch of the twist 
structure is very long, this phase would appear to be a bi- 
axial smectic. The recently identified biaxial SmA phase 
in mixtures of banana and rod-like mesogens may 
correspond to the very long-pitch Sm-AJj (Nt + 2) phase. 
An alternative version of a smectic (Nt + 2)* phase, the 
Sm J 4* L (A r r + 2) phase has n in the plane of the layers as 
depicted for a single layer in Fig. |3^(c) . Since this phase 
is chiral, n will rotate in a helical fashion from layer to 
layer. A more complex (Nt + 2) smectic-C like phase, 
with n making an angle other than or tt/2 with respect 
to N are also possible. 



VII. SUMMARY AND CONCLUSIONS 

In this paper, we have presented a comprehensive 
study of liquid phases of achiral bent-core liquid-crystal 
molecules. Using symmetry we enumerated all possible 
orientationally-ordered liquid phases, classified them by 
subgroups of the rotation group 0(3) under which they 
are invariant, and constructed Landau mean-field theory 
describing these phases and transitions between them. 
One primary conclusion of our work is that in addition 
to the vector (p 1 ) and second-rank nematic (Q 1 -*) order 
parameters, a third-rank tensor order parameter T l ^ k , 
representing third-mass moment, is necessary to char- 
acterize novel liquid-crystal phases of banana molecules, 
such as for example the Nt phase with Did symmetry, 
the tetrahedratic T phase, and the spontaneously chiral 
nematic (Nt+2)* and its chiral polar analog (Vr+2)* . In 
these phases the chiral symmetry is spontaneously bro- 
ken by "condensation" of the biaxial B\i and tetrahe- 
dratic Tq^ order parameters with a nonvanishing angle 
< S < n/2 between their respective principle axes m — 1 



(see Fig.|(b)). 

The Nt phase is neither uniaxial nor biaxial but in- 
stead exhibits an invariance with respect to a four-fold 
improper rotation consisting of a rotation through 7r/2 
about the z-axis followed by an reflection z — ► — z. The 
T phase is invariant under the Td symmetry group of a 
tetrahedron. Like the isotropic phase, its second rank di- 
electric tensor is isotropic, making it optically isotropic, 
but unlike the isotropic phase, it has a nonvanishing sec- 
ond order nonlinear susceptibility Xijl ~ T l ^ k such that 
there is a second order contribution to the polarization 
P^ = Xi%EjEk where Ei is the electric field. 

The chiral nematic (Nt + 2)* phase, like the tradi- 
tional chiral cholesteric and blue phases, will exhibit pe- 
riodic spatial modulations of the direction of molecular 
alignment. Unlike the transition from the isotropic to 
the cholesteric phase, the transitions to the (Nt + 2)* 
phase (from the biaxial N + 2 phase, the tetrahedratic T 
phase and the Nt phase with Did symmetry) are second- 
order, at least in mean-field theory. The pitch of the 
cholesteric structure of the (Nt + 2)* phase diverges as 
these second-order transitions are approached and thus 
changes rapidly with temperature. Since the (Nt + 2)* 
phase spontaneously breaks chiral symmetry, the state 
initially formed upon cooling from the higher-symmetry 
phase will consist of domains of opposite chirality sepa- 
rated by domain wall that will coarsen over time. Chiral 
dopants (or distortions of bent-core mesogens to make 
them chiral) render all phases chiral and, in particular, 
induce a cholesteric pitch of a particular sign on the chi- 
ral extensions (N + 2)*, T* , and N T , of the (N + 2), T, 
and Nt phases. Thus chiral dopants act like an external 
field in an Ising ferromagnet, favoring a particular sign 
of chirality (rather than a particular sign of spin), and 
the transitions from the (TV + 2)*, T*, and N T phases 
to the (Nt + 2)* phase will be analogous to the Ising 
transition in an external magnetic field. In principle, for 
sufficiently large chirality, blue phases, with two orthog- 
onal twist axes can also appear in (Nt + 2)*. In addi- 
tion to these properties of the nonpolar (Nt + 2)* phase, 
the chiral polar (Vt + 2)* phase will exhibit spontaneous 
ferroelectricity, a liquid state that has been a holy-grail 
in liquid-crystal research dating back to Louis Pasteur. 
Light-scattering, circular dichroism and switching with a 
weak electric field would be natural experimental probes 
for these spontaneously chiral states. 

These chiral phases are particularly interesting because 
their smectic analogs, the chiral ShiCU-Pa, SmC s Pa, 
SmC A P F , and SmC s P F (the four B 2 phases) fUJi have 
been realized in banana liquid crystal molecules, gen- 
erating significant excitement in the ferroelectric liquid- 
crystal community. Our work suggests that smectic po- 
sitional order is not necessary and that spatially homo- 
geneous spontaneously chiral liquid-crystal phases are 
generically possible. 

Although, most JT~6[ | experimental systems of banana 
molecules studied so far, appear to undergo direct first- 
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order transitions into smectic phases, our work suggests 
that this situation does not have to be the case, and a 
rich phase structure and hierarchy of continuous transi- 
tions studied here is possible. It is our hope that results 
presented here will stimulate searches in experiments 
and simulations for banana materials that exhibit novel 
orientationally-ordered liquid phases predicted here. 



ACKNOWLEDGMENTS 

Leo Radzihovsky was supported by the National Sci- 
ence Foundation through MRSEC Program at the Uni- 
versity of Colorado at Boulder under #DMR-9809555, 
and by the A. P. Sloan and David and Lucille Packard 
Foundations, and acknowledges hospitality of Harvard's 
Department of Physics, where this work was completed. 
Tom Lubensky was supported by the NSF through grant 
DMR00-96531. The authors thank Noel Clark and Dar- 
ren Link for numerous discussions. 



[1] R.B. Meyer, L. Liebert, L. Strzelecki, P. Keller, J. Phys. 

Lett. (Paris) 36, 69 (1975). 
[2] T. Niori, T. Sekine, J. Watanabe, T. Furukawa, and H. 

Takezoe, J. Mater. Chem. 6, 1231 (1996). 
[3] D. R. Link, G. Natale, R. Shao, J. E. Maclennan, N.A. 

Clark, E. Korblova, D. M. Walba, Science 278, 1924 

(1997). 

[4] In a more faithful nomenclature put forth by the Boul- 
der group in Ref. Q, the so-called B2 phase actually 
refers to four different smectic-C phases SmCiPj of ba- 
nana molecules. These are distinguished by the values of 
the i, j labels. The i label is either S for Synclinic or 
A for Anticlinic indicating, respectively, that molecules 
in neighboring layers tilt in the same or opposite direc- 
tions. The j label is either F for Ferro or A for Anti- 
ferro indicating, respectively that dipoles in neighbor- 
ing layers are either parallel or anti-parallel. The corre- 
sponding four phases are: the interlayer synclinic (same 
tilt from layer to layer) and ferroelectric ShiCsPf (ho- 
mochiral) phase, the interlayer anticlinic and ferroelectric 
SitiCaPf phase, the interlayer synclinic and antiferro- 
electric SvciCsPa phase, and the interlayer anticlinic and 
antiferroelectric StuCaPa phase, all of which have now 
been observed. 

[5] G. Pelzl, S. Diele, W. Weissflog, Adv. Mater. 11, 707 
(1999). 

[6] J. P. Bedel, J. C. Rouillon, J. P. Marcerou, M. Laguerre, 
H. T. Nguyen, and M. F. Achard, Liq. Cryst. 28, 1285 
(2001). 

[7] D. M. Walba, E. Korblova, R. Shao, J. E. Maclennan, 
D. R. Link, M. A. Glaser, and N. A. Clark, Science 288, 
2181 (2000). 

[8] Banana-Shaped Liquid Crystals: Chirality by Achiral 



Molecules, Workshop, Technical University of Berlin, Afb 

335, Anisotrope Fluid' (1997). 
[9] H. R. Brand, P. E. Cladis, and H. Pleiner, Europhys. J. 

B 6, 347 (1998). 
[10] A. Roy, N. V. Madhusudana, P. Toledano, A. M. 

Figueiredo Neto, Phys. Rev. Lett. 82, 1466 (1999). 
[11] M. P. Neal, A. P. Parker, and C. M. Care, Mol. Phys. 

91, 603 (1997). 

[12] P. J. Camp, M. P. Allen, and A. J. Masters, J. Chem. 

Phys. Ill, 9871 (1999). 
[13] J. L. Billeter and R. A. Pelcovits, Liq. Cryst. 27, 1151 

(2000). 

[14] R. Memmer, Liq. Cryst. 29, 483 (2002). 

[15] Y. Lansac, P. K. M aiti, N. A. Clark, and M. A. Glaser, 
|cond-mat/0202154 . 

[16] K. Kishikawa, et. al. Chem. Matt. 11, 867 (1999); D. 
Shen, s. Diele, G. Pelzl, I. Wirth, and C. Tschierske, J. 
Mater. Chem. 9, 661 (1999); J. Matraszek, et al, Liq. 
Cryst. 27 429 (2000); 

[17] T. C. Lubensky and L. Radzihovsky, unpublished. 

[18] See for example Michael Tinkham, Group Theory and 
Quantum Mechanics (McGRaw-Hill, N.Y., 1964) 

[19] John Toner, Phys. Rev. A 27, 1157-63 (1983) 

[20] By "anisotropic" we mean an orientationally-ordered 
phase, which, as in a case of the tetrahedratic phase 
T, characterized by S — and isotropic dielectric ten- 
sor, should appear isotropic in linear optics experiments. 
It is nevertheless an anisotropic phase since its rota- 
tional symmetry is smaller than that of fully disordered 
isotropic liquid, a distinction that can in principle be ex- 
perimentally observed in nonlinear optical experiments, 
such as optical second harmonic generation. 

[21] P. G. de Gennes and J. Prost, "The Physics of Liquid 
Crystals", Clarendon Press, Oxford (1993). 

[22] P. Chaikin and T. C. Lubensky, "Principles of Con- 
densed Matter Physics," Cambridge University Press, 
Cambridge (1995). 

[23] R. Alben, Phys. Rev. Lett. 30, 778-81 (1973). 

[24] In mean-field theory, a noncritical order parameter ex- 
plicitly induced by spontaneous ordering of another order 
parameter has its growth with reduced temperature and 
field, completely determined by the critical exponents 
controlling the behavior of the critical order parameter. 
Canonical example is an XY-model where ordering of 
V>i explicitly induces higher harmonic order parameters, 
such as e.g., ip n through coupling ipiip^Thechiral + c.c. 
However, despite of the irrelevance of this coupling for 
the critical behavior of ^1 (noncritical can be simply 
integrated out, finitely renormalizing bare parameters of 
the H[tpi]), because of the nonlinear nature of this cou- 
pling operator, beyond mean-field theory, it, as well as 
susceptibility \n for ip n generically receive nontrivial fluc- 
tuation corrections and therefore have independent and 
nontrivial scaling dimensions. These can be shown to lead 
pB] to breakdown of mean-field f3 n = n/3i relation, as is 
most strikingly seen for an XY-model in two dimensions. 

[25] A. Aharony, R.J. Birgeneau, J.D. Brock, and J.D. Lit- 
ster, Phys. Rev. Lett. 57, 1012 (1986) 

[26] LG. Fel, Phys. Rev. E 702, 702 (1995). 

[27] In the preceding reference, Fel studied the I — > T tran- 



27 



[28] 
[29] 

[30] 



[31] 
[32] 
[33] 



sitions and many of the interesting properties of the T 
phase using a Landau theory for a third-rank symmetric- 
traceless tensor order parameter. He failed, however, to 
identify all of the invariants in the Landau theory and, 
therefore, failed to identify the / — » N + 3 transition. He 
also missed the second quartic invariant in r lJ (the vt 
term in our Eq.4.4), and therefore incorrectly identified 
the I — > T transition as being in the 0(7) universality 
class. For correct analysis, see Ref. [^ 
R.G. Priest and T.C. Lubensky, Phys. Rev. B 13, 4159 
(1976). 

Thermal fluctuations drive the I-T transition first-order 
as shown in L. Radzihovsky and T. C. Lubensky, Euro. 
Phys. Lett. 54, 206 (2001). 

The reduced transition temperatures for ordering from 
the nematic into polar states with one of the three order 
parameters pz, p2, and p\ finite are respectively deter- 



mined by r, 



r eff 



r p + 4w p qS/3, 



eff 



2w p qS/3, 



•~ p — 2w p qS/3. Nematic ordering along n, there- 
fore obviously raises T c for pz, and lowers it equally for 
pi and pi. 

K. G. Wilson and J. Kogut, Phys. Rep. C, 12, 77 (1977). 
J. Rudnick, Phys. Rev. B, 18, 1406 (1978). 
The nature of the ordered phase is determined by the 
sign of the coefficient of the cubic symmetry breaking 
term. If positive, the XY-order parameter pi t 2 will order 
along either the m or 1 axis, leading to pi t 2 = pi(l,0) 
or p\^ = P2(0, 1) order. For negative value of the cubic 
anisotropy coefficient, the ordering is along the diagonal 



[34] 



[35] 
[36] 



with pi.2 — p(l, 1). With the Nt phase characterized by 
nonzero Tj and Ta — 0, it is only the latter (1,1) ordering 
that leads to the N + V phase. In the former case of (1,0) 
or (0, 1) order, the resulting state is the chiral (N + V)* 
state characterized by triclinic Ci symmetry. Evidence 
for the smectic version of this state, the so-called SmCc , 
proposed by de Gennes |2lj has recently been presented 
in freely suspended films and in the bulk, and identified 
with the £?7 texture. |m| 

D. M. Walba, E. Korblova, R. Shao, D. N. Coleman, 
N. Chattham, and N. A. Clark, in Proceedings of the 
ALCOM Symposium, Cuyahoga Falls, Ohio, (2000), Ab- 
stracts, Vol. XII, and the Eighth International Confer- 
ence on Ferroelectric Liquid Crystals (2001), 0-2; A. 
Jakli, D. Kruerke, H. Sawade, and G. Heppke, Phys. Rev. 
Lett, 86, 5715 (2001). 

E. Brezin, LeGuillou and J. Zinn- Justin, Phys. Rev. B, 
10, 892 (1974). 

It is easy to show using Eq. (|5.22h for f^ JV) that within 



[37] 



the quartic /q model, Eq. (4.3), valid only close to the 
I-N transition (and even then, only if it is weakly first- 
order, so that S is small), the — 6wqS is strictly 
positive, precluding an iV — > N + 2 transition. However, 
a more general free energy density of Q lJ (e.g., includ- 
ing higher order nonlinearities in Q lJ ) can be shown to 
lead to fg^ that does go negative, thereby triggering 
the transition to a biaxial N + 2 state. p!j| , pl| 
R. Pratibha, N. V. Madhusudana, B. K. Sadashiva, Sci- 
ence 288, 2184 (2000). 



28 



